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ABSTRACT 


The  theory  and  assumptions  used  in  developing  equations  for  the  error 
analysis  of  a  general  class  of  inertial  navigation  systems  are  described. 
The  computer  program  developed  for  their  solution  is  described  from  a 
user's  point  of  view.  Its  application  includes  the  synthesis  and/or 
analysis  of  inertial  navigation  systems  used  in  ballistic  missile  or 
terrestrial  space  missions.  The  program  is  designed  to  allow  studies 
of  both  pure  inertial  and  aided  inertial  navigation  systems,  the  latter 
being  the  process  of  updating  navigation  data  via  data  from  external 


sensors. 
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SECTION  1 


INTRODUCTION 


This  report  describes,  from  a  user's  point  of  view,  a  computer  program 
developed  for  the  error  analysis  of  a  general  class  of  inertial  navigation 
systems.  It  is  generally  understood  that  these  systems  will  be  used  in 
connection  with  either  ballistic  missile  or  space  mission  applications. 

Section  2  describes  the  classes  of  inertial  navigation  systems  considered, 
develops  the  equations  necessary  to  perform  an  error  analysis,  enumerates 
the  assumptions  made  in  their  derivations,  and  describes  the  methods  and 
equations  used  for  data  presentation.  The  equations  in  Section  2  form  the 
bases  for  the  computer  program,  which  was  developed  to  solve  them. 

Section  3  deals  with  the  operational  aspects  of  using  the  computer  program 
to  perform  error  analyses.  The  input  data  requirements  and  procedures 
are  discussed  and  the  resulting  output  data  and  formats  described.  The 
logical  order  of  the  computations  resulted  in  the  development  of  two 
independent  programs:  The  first,  called  ERAN,  solves  the  equations 
presented  in  Section  2.  3;  the  second,  called  0 UTP,  solves  those  presented 
in  Section  2.  4.  These  were  progr  mmed.  for  the  IBM  7090/7094  to  be  rim 
under  the  control  of  the  IBM  Basic  Monitor  (IBSYS)  Programming  System 
with  the  assumption  that  core  is  set  to  zero  before  loading  of  the  programs. 
There  is  a  certain  amount  of  intentional  redundancy  in  the  material  pre¬ 
sented  in  Sections  2  and  3.  This  was  done  so  that  once  one  ir,  famiiia~  with 
the  material  presented  in  Section  2,  it  will  only  be  necessary  to  refer  to 
Section  3  for  program  operations. 

Section  4  presents  three  sample  test  cases  illustrating  the  input  data  pro¬ 
cedures  and  the  formats  of  output  data,  and  demonstrating  some  of  the 
flexibility  and  capabilities  of  the  program. 


Appendices  A  through  C  present  material  augmenting  the  mam  body  of  the 
report.  Appendix  A  consists  of  Standard  Input  Data  Forms,  Appendix  B 
contains  EK AN  and  0UTP  input  Data  for  Sample  Cases,  and  Appencjjix  C 
gives  the  Output  Listings  for  Sample  Cases. 


Appendix  D  is  devoted  to  the  subject  of  updating  or  correcting  navigation 
data  through  the  use  of  external  data  sources  utilizing  various  sensor 
configurations.  Algorithms  are  derived  for  three  possible  schemes  of 
data  processing.  Although  equations  have  not  been  programmed,  the 
logical  structure  of  the  computer  program  is  designed  so  that  this  fea^re 
can  be  incorporated  without  major  revisions. 


/ 


In  Appendix  E  the  method  is  discussed  of  treating  the  effects  of  aero¬ 
dynamic  drag  in  orbit  when  the  accelerometers  are  disconnected  from 
the  navigation  computer. 


Appendix  F  contains  all  the  figures  called  out  in  the  report  and  Appendix  G 
the  program  definitions  and  constants. 


SECTION  2 


EQUATION  DEVELOPMENT 

2.  1  INTRODUCTION 

This  section  defines  the  classes  of  inertial  navigation  systems  considered 
and  develops  the  equations  necessary  to  perform  an  error  analysis  of  a 
given  configuration. 

Section  2.  3  relates  the  derivation  of  the  differential  equations  of  navigation 
error  to  a  broad  class  of  system  errors.  The  solution  of  these  equations 
results  in  the  linear  transformations  (sensitivities)  of  navigator  errors 
into  errors  of  navigation  data.  The  classes  of  errors  include  those  of 
sensor  anomalies,  initial  conditions,  terminal  control,  and,  for  certain 
operational  procedures,  the  effect  of  orbital  drag  uncertainties. 

Section  2.4  describes  the  equations  used  for  processing  these  sensitivities 
into  individual  navigation  vector  errors  for  each  error  source,  and  those 
which  statistically  sum  all  vector  errors  presented  as  a  covariance  matrix. 
Various  coordinate  systems  and  processing  methods  are  described  for  data 
presentation. 
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2.2  NAVIGATION  SYSTEM  CONFIGURATIONS 

Th«  inertial  navigation  system  configurations  considered  are  schematically 
presented  in  Figure  1.  The  equations  developed  for  error  analysis  are 
applicable  to  torqued  or  inertially  oriented  gimballed- platform  systems, 
and  to  a  certain  class  of  strapped-down  systems.  The  essential  sensors 
used  by  the  navigation  system  are  three  accelerometers,  which  sense  the 
magnitude  of  the  applied  external  accelerations,  and  three  gyros,*  which 
?<*nse  angular  dynamics  so  that  the  direction  of  the  applied  acceleration 
can  be  derived.  The  constraints  of  accelerometer  mounting  are  such  that 
the  three  sensing  (input)  axes  cannot  be  coplanar,  but  can  be  nonorthogonal. 
Gyros  are  assumed  to  be  mounted  in  a  triad  so  that  their  sensing  (input) 
axes  are  orthogonal.  The  method  of  derivirg  accelerometer  orientation 
from  gyro  signals  is  assumed  to  be  one  of  the  following  three: 

a.  Gimballed  Platforms.  In  this  configuration  -  the  most 
conventional  -  the  platform  is  initially  aligned  to  some 
auxiliary  references.  For  initial  alignment  on  the  ground, 
the  gravity  vector,  which  is  sensed  by  pendulums  or  the 
accelerometers,  is  used  for  vertical  reference;  azimuth 
is  referenced  to  either  optical  sensors  or  a  gyro  com¬ 
pass.  For  in-orbit  alignment,  stellar  or  horizon  sensors 
are  used.  The  gyros  measure  any  deviation  of  the  plat¬ 
form  from  the  initial  alignment  and  their  signals  are 
used  to  torque  the  platform  in  such  a  way  that  they  be¬ 
come  null,  thus  maintaining  the  initial  reference.  In 
some  cases,  the  platform  is  torqued  either  to  reduce  the 
total  gimbal-  angle  travel,  or  to  maintain  prelaunch 
(earth)  rates.  To  achieve  this,  the  gyros  are  torqued  at 
prescribed  rates,  which  the  platform  follows.  In  this 
case,  the  accelerometer  transformation  matrix  is  a 
function  of  time  and  computed  from  the  commanded  rates. 
Mathematically,  this  is  identical  to  configuration  (b). 

b.  Strapped- down/ Caged  Gyros.  In  this  configuration,  the 
platform  is  mounted  directly  (possibly  by  a  shock  mount) 
to  the  airframe.  The  gyros  sense  angular  deviations 
from  the  initial  reference,  but  are  torqued  to  null  their 


* 


*Two  single-degree-of-freedom  gyros  can  be  oriented  to  represent  one 
two- degree- of- freedom  gyro. 
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signals.  These  torquing  signals  are  a  measure  of  the 
rate- of- change  of  the  platform  orientation  and  are  used 
to  compute  the  transformation  matrix  required  for  the 
accelerometers.  The  algorithm  used  is  based  on  the 
matrix  differential  equation  of  direction  cosines. 

c.  Strapped- down/ Free  Gyros.  Until  the  adyent  of  the 
electrostatic  gyro  (fiSu),  lighweight  wide-  angle  free 
gyros  were  not  sufficiently  accurate  to  be  considered 
for  this  application.  Free  (two- degree-of- freedom) 
gyros  are  used  in  platform  configurations  but  are  re¬ 
stricted  to  small- angle  deviations.  With  its  high  degree 
of  accuracy,  the  ESG  acts  as  a  potential  attitude  reference 
sensor.  In  this  configuration,  the  angles  the  gyro  case 
(thus,  the  accelerometer)  makes  with  respect  "o  the  spin 
axis  of  the  gyro  are  read  out  and  used  to  compute  the 
transformation  matrix. 

The  sensor  data  is  processed  by  a  computer  to  derive  position  and  velocity 
data.  It  is  assumed  that  the  navigation  system  computer  has  perfect 
algorithms  for  gravity,  and  that  its  word  length  and  integration  schemes 
are  such  as  to  produce  negligible  errors.  Most  accelerometer  outputs  are 
in  the  form  of  pulse  rates  proportional  to  acceleration,  thus  there  are 
additional  complications  in  deriving  inertial  velocity  when  the  platform  is 
not  inertially  oriented.  The  algorithms  used  in  these  cases  become  the 
subject  of  a  separate  analysis,  which  requires  detailed  knowledge  of  the 
hardware  characteristics.  In  general,  a  special- j  irpose  computer  (e.  g, , 

,a  Digital  Differential  Analyzer)  would' be  required  to  buffer  the  processing 

« 

of  accelerometer  and  gyro  data  into  a  suitable  form  for  processing  by  a 
general-purpose  computer.  Generally,  this  form  is  sensed  velocity  data, 
which  is  then  corrected  for  the  effects  of  gravity  from  which  the  trajectory 
position  and  velocity  data  are  determined.  In  this  analysis  it  is  assumed 
that  the  error  in  these  computations  small  enough  to  be  negligible,  or 
convertible  to  equivalent  sensor  errors;  therefore,  the  errors  in  indicated 
position  and  velocity  are  functions  of  sensor  anomalies  only. 

Provisions  are  included  for  analyses  of  aided  inertial  navigation  systems 
in  which  external  sensors  are  used  for  measuring  position,  velocity,  and/or 


platform  orientation.  The  measurements  are  processed  by  various  techni¬ 
ques  and  applied  as  corrections  to  the  navigation  system  data.  These 
corrections  are  discussed  in  Appendix  D  under  Resets.  The  computer 
program  in  its  present  form  does  not  include  the  coding  of  these  equations. 


2.3 


ERROR  SENSITIVITY  EQUATION  DEVELOPMENT 


2.3.1  Coordinate  Systems  and  Transformation  Matrices 

The  basic  coordinate  system  used  for  computing  navigation  errors  is  an 
earth-centered  inertial  (ECI)  system,  in  which  the  Z  axis  is  along  the  earth's 
polar  axis,  and  the  X  and  Y  axes  lie  in  the  earth's  equatorial  plane,  forming 
a  right-hand  orthogonal  axis  system.  Generally,  the  convention  is  that  the 
X  axis  passes  through  the  Greenwich  meridian  at  time  zero. 


The  notation  used  for  coordinate  transformation  matrices  is  the  symbol  M 
with  two -lettered  subscripts  to  identify  the  respective  coordinate  systems; 
e.g. ,  identifies  the  transformation  matrix,  which  is  used  to  transform 

vectors  to  the  ECI  coordinate  system  from  the  K  coordinate  system.  Con¬ 
versely,  transforms  vectors  to  the  K  system  from  the  ECI  system.  The 

necessary  coordinate  systems  and  the  transformation  matrices  are  described 
in  the  following  paragraphs. 


2.3. 1.1  Platform  Coordinate  System 

The  platform  axes  are  designated  Pj,  P^,  and  P^  and  form  a  right-hand 
orthogonal  coordinate  system.  The  initial  transformation  matrix  is  developed 
by  considering  the  platform  axes  to  be  initially  aligned  with  the  ECI  axes  and 
by  applying  ordered  rotations  of  0p  about  P^,  a  negative  \p  about  P^,  and  a 
negative  4^p  about  P^.  Thus,  the  matrix  that  transforms  vectors  in  ECI  co¬ 
ordinates  to  vectors  in  platform  coordinates  is 
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Figure  2  illustrates  the  usual  initial  platform  orientation,  where  0p  is  the 
longitude,  \p  is  the  geocentric  latitude  (positive  North  latitude),  and  ‘J'p  is 
the  aainmth  (positive  is  conventional  asimuth  from  North). 

Sinco  pUtfcrm  coordinates  are  orthogonal 


mep  =  mpe 


u 


t 


where  T  denotes  the  transpose.  Mpg  can  be  a  function  of  time  (for  strapped  - 
down  or  torqued  platforms)  and  its  calculation  is  discussed  in  Section  2.3.2. 

2.3. 1.2  Gyro  Coordinate  System 

It  is  necessary  to  assign  a  coordinate  system  to  each  gyro  so  that  the  gyro 
errors  can  be  determined.  The  gyro  axes  for  each  component  are  right- 
hand  orthogonals  and  designated  °il*  Ci2  ,  and  G.j  (i  =  number  1,  2,  or 
3  gyro),  where  Gjj  is  the  sensing  (input)  axis  of  the  ith  gyro.  It  is  assumed 
that  Gjj,  G21,  and  G^j  also  form  a  right-hand  orthogonal  axis  system. 

As  a  result  of  this  assumption,  the  development  of  the  gyro  coordinate  sys¬ 
tems  is  simplified.  Since,  for  any  one  configuration,  the  gyro  axes  are 
assumed  to  be  fixed  with  respect  to  the  platform  axes,  the  gyro  coordinate 
systems  are  developed  with  respect  to  the  platform  coordinate  system.  Fig¬ 
ure  3  illustrates  the  initial  gyro  alignments  with  respect  to  the  platform  axes. 
The  method  of  specifying  gyro  orientations  is  by  specification  of  an  axis  (1, 

2,  or  3)  and  an  argument  (angle*)  of  successive  rotations.  Each  rotation 
operates  on  the  gyros  as  a  triad;  i.e. ,  all  three  gyros  are  being  rotated  and 
thus  are  maintaining  their  axis  c.-ientation  with  respect  to  each  other  during 
the  rotations.  The  axis  of  rotation  referred  to  above  is  that  of  the  number 
one  gyro.  Thus,  the  matrix  that  transforms  vectors  in  platform  coordinates 
to  gyro  coordinates  (gyro  input  axes)  is 


M, 


GP/ 
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( 


where 


and 


i  =  the  axis  of  rotation  (i  =1,  2,  or  3) 
th 

j  =  the  j  rotation  (j  =  lj2,  ...  up  to  5) 

a-*-. 


Tir 


T 

2j 


‘ce. 


.  se- 

L  J 


ce. 

se. 

j 

j 

-se. 

ce. 

j 

j 

0 

-se. 

j 

l 

0 

0 

ce. 

JJ 


transformation  for  a 
rotation  about  axis 


transformation  for  a 
rotation  about  axis 


se. 

J 

ce. 

J 


0* 

0 


0 


transformation  for  a 
rotation  about  G^  axis 


where  0.  is  the  angle  of  the  j^1  rotation  and  a  positive  angle  is  in  the  sense  of  a 
right-hand  rotation.  Upon  completion  of  this  set  of  rotations,  there  remains 
an  additional  degree  of  rotational  freedom  of  each  gyro  about  its  input  axis. 

By  considering  this  degree  of  freedom,  the  vector  components  along  the  other 
two  axes  of  the  gyro  (axes  2  and  3)  are  determined,  utilizing  the  matrices 


'  0 

S4-J* 

MG2  = 

S*2 

0 

C*2 

-c*3 

S*3 

o  . 

transforms  a  vector  in 
gyro  coordinates  to  vector 
components  along  the  2  axis 
of  each  gyro 


'  0 

-S4»j 

CvJ/j" 

MG3  = 

c*. 

0 

■S4*2 

transforms  a  vector  in  gyro 
coordinates  to  vector  com¬ 

£• 

ponents  along  the  3  axis  of 

.’54*3 

C»)<3 

0 

each  component 

where  is  the  angle  of  rotation  for  the  i^1  gyro  and  a  positive  angle  is  in 
the  «ense  of  a  right-hand  rotation. 

It  is  convenient  to  describe  the  gyro  axes  in  the  model  of  a  single-degree -of - 

freedom  gyro,  where  the  G.  j  is  the  input  reference  axis,  is  the  output 

axis,  and  G.,  is  the  spin  reference  axis.  Then  a  rotation  of  y.  =  90  deg  of  ' 
th  i 

1  gyro  will  provide  an  orientation  of  two  gyros  with  orthogonal  input  axes 

and  parallel  spin  axes,  a  model  of  a  two -degree -of-freedom  gyro.  Since  the 
gyro  (input  axes)  coordinate  system  is  orthogonal,  the  following  transforma¬ 
tions  are  derived 


MPG  =  mgp 
mge  =  mgpmpe 

meg  =  ^mgpmpe^  =  mepmpg 

2.3.  1.3  Accelerometer  Coordinate  System 

The  following  two  options  are  used  for  specifying  the  alignment  of 
accelerometers . 

FIRST  OPTION 

First  is  the  specification  of  an  orthogonal  triad  and  the  method  of  specifica¬ 
tion  is  identical  with  that  described  for  the  gyro  components;  that  is,  an  axis 
and  argument  of  successive  rotations  are  specified  to  align  the  accelerometer 
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SECOND  OPTION 


The  second  option  allows  for  u  non  orthogonal  accelerometer  configuration. 
Here  the  method  of  specification  iw  the  same  (i.e. ,  specification  of  an  axle 
and  an  argument);  however,  each  accelerometer  is  specified  independently. 
The  initial  orientation  cf  each  accelerometer  is  the  same  with  the  accelerom¬ 
eter's  input  Aiiis  aligned  with  Pj,  with  P^,  and  A^  with  P^.  There¬ 
fore,  a  matrix  is  developed  for  each  accelerometer  k,  as  follows: 


M 


APk 


=  T. 


ijk 


T:'.2kTilk  ^  =  ** 


2,  and  3) 


From  these  three  matrices,  the  matrices  M^p,  M^,  and  are  formed 


as 


AP 


M 


M 


11 


12 


UMn 


M 


A?. 


M 


21 


M 


22 


LM23j 


M 


PA 


MA3  = 


M 


31 


1 


Mr 


LM33^ 


M 


PA 


where  M.,  is  the  i  n  row  of  M^p^  (k  =  1,  2,  3) 


"ik 
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Since  for  this  option  the  accelerometers  are  nonorthogonal,  Mp^  i*  developed 
from  the  inverse  rather  than  the  transpose,  and  the  following  relationships 
result 


MPA  =  MAP 


-1 


mae  =  mapmpe 


MEA  ^“ApMpE1'1  ■  MEFMPA 
2.3. 1.4  Initial  Condition  Coordinate  System 

These  coordinate  systems  are  used  to  transform  errors  specified  in  a  con¬ 
venient  coordinate  system  into  errors  in  ECI  coordinates  .  The  initial  platform 
errors  are  specified  in  the  platform  coordinate  system  and  transformed  into 
ECI  coordinates  by  Mgp-  Two  options  are  provided  for  specifying  initial 
position  and  velocity  errors.  In  the  first,  the  errors  are  assumed  to  be 
referenced  to  platform  axes;  thus  M£p  is  used.  In  the  second,  the  errors  are 
assumed  to  be  referenced  to  the  geocentric  vertical  and  directed  in  azimuth 
by  a  specified  angle  (4>j)  Figure  5).  Tbs  transformation  matrices  for 
position  and  vp1ocity  errors  are  then 

OPTION  1 


^Ei  =  mep 
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OPTION  2 


C0  -S0 
90  C0 


CX  0  -SX 

0  1  0 


SX  0  CX 


0 

0 

C*! 

S4*j 

-S*J 

where 


S0-^ 


C0  = 


-fi 


Vy2 


CX=| 


R  =  JD2  +  Z* 


2.3. 1.5 


X,  Y,  Z  -  initial  (t  =  tQ)  ECI  position  coordinates 
0  and  X  =  the  initial  longitude  and  geocentric  latitude,  respectively 
4^  =  the  azimuth  orientation  of  the  coordinate  system. 
Terminal- condition  Coordinate  System 


These  coordinate  systems  are  used  when  it  is  desired  to  propagate  the  effects 
of  terminal-control  errors  generated  by  other  systems  for  which  there  can  be 
no  further  corrections  (e.g.,  guidance  steering  and  thrust  tailoff  errors  at 
thrust  termination  of  a  ballistic  missile).  The  transformation  for  terminal - 
position  errors  is 


M 


ET1 


C0  -S0 
50  C0 
0  0 


0 

0 

1J 


C>  0  -sx 

0  1  0 

Lsx  o  cx. 


1  0  0 

0  Cv|i  S4< 

Lo  -S4»  c*J 


where 


SX  = 


S4>  = 


*5 


cx=| 


VN 

C4<  =  rr- 

VH 


D=J 


X2  +  Y2 


R  =s/D27z2 


«  • 


„  _  -YX  +  XY 

VE  - - B - 


„  _  -Z(XX  +  YY)+  ZD* 

VN  = - KB - - 


H 


=  /v2  +  V2 

n/  N  VE 


and 

X,  Y,  Z  =  terminal  ECI  position  coordinates 

X  Y  Z  =  terminal  ECI  velocity  coordinates 

Vg  .=  east.,  north,  and  hor  izontal  velocity 

vector  components 

0,  X,  <j>  =  longitude,  latitude,  ar.d  velocity  vector 
azimuth  at  the  termination  time 
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The  transformation  for  terminal  velocity  errors  is 


where 


'  Cy 

0 

Sy* 

**1X2  *  MET1 

0 

1 

0 

.-Sy 

0 

Cy. 

„  XX  +  YY  +  ZZ 

SY  - - RV - 


V  =  (X2  +  Y2  +  Z2)112 


and 

y  is  the  flight  path  angle  and  equals  the  angle  that  the 
velocity  vector  makes  with  the  horizontal  plane 

V  is  the  magnitude  of  the  velocity  vector. 

These  coordinate  systems  are  illustrated  in  Figure  6. 

The  other  coordinate  systems,  used  for  output  data  and  reset,  are  described 
in  Section  2.  4.  1  and  Appendix  D. 

2.  3.  2  Navigation  Kinematics 

For  purposes  of  error  analyses,  it  is  of  interest  to  derive  the-  sensitivities  of 
navigation  and  platform  orientation  errors  as  a.  function  of  the  anomalies  of 
system  parameters  (which  are  random  variables).  These  sensitivities  are 
then  used  with  the  statistics  of  the  system  parameters  to  develop  the  statistical 
characteristics  of  the  navigation  data.  For  this  purpose,  the  differential  equa¬ 
tions  are  linearized  so  that  all  the  advantages  of  linear  analysis  can  be  utilized. 
The  following  sections  develop  this  procedure. 


2.  3.  2.  1  Equations  of  Motion 

The  differential  equation,  expressed  in  ECI  coordinates,  that  describes  the 


vehicle  equations  of  motion  for  an  assumed  point  mass  system  is 


RJ  * 


where 


p  =  gravitational  constant 

R  =  magnitude -of-position  vector 

G  =  vector  of  higher -order  gravity  terms 
•  * 

X  =  acceleration  vector  of  the  vehicle 


Xg  =  sensed  acceleration  vector  due  to  all  external  forces 
(thrust,  drag,  etc.  ) 


The  navigation  system's  computer  mechanizes  and  solves  this  equation  with 

*  ♦ 

the  appropriate  initial  conditions  of  position  and  velocity  (Xq  and  Xq).  The 


sensed  acceleration  components  are  measured  by  the  accelerometers  and 
transformed  into  ECI  coordinates  resulting  in 


X  =(MpoMd.],A 
sm  1  EP  PAJ 1  i 


oc 


where 

♦  • 

X  =  the  sensed  acceleration  measured  and 
sm 

transformed  into  ECI  coordinates 

[Mj,pMpA]i  =  the  transformation  between  accelerometer  and 
ECI  coordinates,  which  the  computer  uses 

A  =  accelerometer  measurements  corrected  for 
calibrated  anomalies 


The  choice  of  coordinate  systems  made  here  was  for  convenience.  The 
results  are  invariant  with  the  computational  coordinate  system  chosen, 
provided  that  the  assumption  of  perfect  computer  algorithms  is  valid. 


The  corrected  accelerometer  measurements  are  the  result  of  input  acceler 
ations  and  functionally  related  correction  terms.  Thus,  the  accelerometer 
equation  becomes 


Xoc  =*,+  >■(*> 

«iMApMpE]2X<t/<X) 


where 


Ag  =  uncorrected  accelerometer  measurements 

/(A)  -  the  equation  for  the  accelerometer  anomalies,  which 
includes  corrections  for  biases,  scale  factor,  non¬ 
linearity,  etc.  ,  based  on  instrument  calibrations. 

(This  function  is  discussed  in  Section  2.  3.  3, 

Error  Sources) 

[M.pMpE]_  =  the  transformation  between  actual  accelerometer  axes  and 
ECI  coordinates,  a  function  of  platform  errors,  acceler¬ 
ometer  alignments,  nonlinearities,  etc. 

Therefore,  the  equation  that  is  solved  by  the  navigation  computer  is 


X  =  -  -^[X  +  G]  +  [MEpMpA]l{[MApMpE]2Xs  +  /(A)} 
R 


Generally,  it  is  not  required  that  a  distinction  be  made  between  [Mep^p^I 
and  [M^p^pE^*  except  in  cases  where  a  given  error  source  affects  both 
computer  and  platform  transformations,  e.  g,  ,  initial  position  errors. 

2.  3.  2.  2  Constraints  and  Assumptions 

Before  proceeding  to  the  linearization  of  the  navigation  system  equations,  the 
following  constraints  are  imposed. 

a.  Nominal  Trajectory  Reference 

Th'*  parameters  that  affect  the  sensed  acceleration  vector 
(e.  g.  ,  thrust,  weight,  wind,  control  system,  etc.  )  are 
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statistically  independent  of  the  navigation  parameters 
(accelerometers,  gyros,  etc.)*  Therefore,  to  include 
them  in  the  analysis  would  result  in  trajectory  deviations 
constrained  by  the  guidance  system,  but  not  necessarily  in 
navigation  errors  of  the  measured  trajectory.  The  only 
way  in  which  the  two  sets  of  parameters  could  be  entered 
into  the  analysis  would  be  from  nonlinearities  for  which 
Montt.  Carlo  techniques  would  be  required  to  develop  their 
effects..  Fortunately,  these  nonlinearities  are  smau  and 
can  be  assessed  independently  by  performing  a  guidance 
error  analysis  and  the  navigation  system  error  analysis, 
separately,  using  the  same  reference  trajectory.  The 
navigation  error  analysis  can  be  repeated  for  the  maid- 
mum  perturbed  trajectory  developed  in  the  guidance  error 
analysis  to  determine  if  further  analysis  is  required. 

b.  First-Order  Partials 

If  the  equations  are  expanded  into  a  Taylor  series,  the 
relative  magnitudes  of  second-order  partials  can  be  deter- 
mined.  For  example,  by  taking  one  component  of  the 
gravity  expression 

X  =  X  =  -  %  (assuming  R  =  X,  Y  =  Z  =  0) 

R  XT 

and  expanding  it  to 


it  is  seen  that  the  maximum  effect  of  the  second- order 
partial  is  a  function  of  AR/R.  Since  AR/R  «<  1  for  any 
reasonable  system,  it  can  be  neglected.  Many  analyses 
of  inertial  guidance  systems  neglect  completely  the  first 
order  gravity  partials*  or  treat  them  as  constants.  ** 


See  for  example  Reference  2,  p  304: 

S' 

See  for  example  Reference  3. 


"Use  of  Normalized  Integrals.  " 
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With  a  •imilftr  analysis,  the  second-order  gravity  terms 
(fi)  can  be  eliminated,  even  though  they  present  linear 
terms  In  the  equations. 

c.  Small-Angle  Approximations 

Applying  the  small-angle  approximation  to  platform  and 
computer  errors  is  justified  on  the  basis  of  comparing 
second-order  partials.  With  this  approximation,  small  - 
angle  rotations  can  be  represented  as  vectors  and  the 
following  vector  matrix  relationships  can  be  used.  * 

(1)  Vector  Transformations 


M. 
ij  J 


where 

4 ^  =  angular  error#  expressed  in  the  i  coordinate 

system  due  to  angular  errors  in  the  j  coordinate 
system  4^ 

M..  =  the  transformation  matrix  to  the  i  coordinate 
J  system  from  the  j  coordinate  system. 

For  example 


r*  1 

X 

\  =  K  = 

4 

y 

mep  = 

1 

•a 

N 

1 _ 

.*3. 

(2)  Matrix  Differentials 


[6M. .].  -  [*.]M.. 
1  lj  1  1  1  ij 


r  0 


■♦3 


*3 

0 

-♦1 


*1 

0  . 


M.. 


See  Reference  4, 


pp  251  to  253. 
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where 


=  the  change  in  Mjj  due  to  small  rotations 
in  the  i  coordinate  system, 

[® .  ]  =  a  skew  symmetric  matrix  formed  from 

the  small  angle  rotations  in  the  i  coordinate 
system. 

For  example 


I  6M..]. 
ij  1 


when  0.  are  small  rotations  in  platform  coordinates. 

Also 


[6MEp]p 


=  f6Mw]£  -M, 


If  the  small  rotations  are  expressed  as  rates  times 
time  (05t  =  w6t),  then  in  the  limit  as  5t  -♦  0,  the 
matrix  differential  equation  is 


f  !  p 

6t 


0 

0 

-w 

1 
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and 


mpe  =  (np]MpE 


( 3 )  Matrix  Transformations 


!$•]  = 

j  ji  iJ  ij 


where 

[$.]  =  skew  symmetric  matrix  expressed  in  the  j 
^  coordinate  system  due  to  rotations  in  the  i 
coordinate  system. 

For  example 

[®E]  =  MEp[$p]MpE 

that  is 


f° 

« 

-0 

y 

r  o 

*3 

-H 

K 

0 

0 

X 

mep 

-0, 

0 

*i 

U 

-p 

X 

0  . 

-  0  ^ 

o  - 

This  can  be  shown  by  using  vector  matrix  relation¬ 
ships  as  follows:  Let  vectors  in  :he  P  coordinate 
system  be  i  .  iated  as 


u 


(. 
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0px  Ap 

■  0 

A3 

1 

> 

tv 

_ i 

=  -A 

PX  ^P 

'  0 

03 

~h 

As 

0 

A1 

*P  =  - 

-03 

0 

+1 

'  A2 

-A1 

0  - 

“  *2 

-01 

0  - 

=  -  [® p]  Ap 

r  r  l 

Transforming  the  above  vectors  into  E  coordinates 
results  in 

=  MEPYP  =  "MEP^P^AP 
*E  =  MEP^P 
ae  =  mepap 

The  vectors  Y,  <j>,  and  A  are  invariant  under  an 
orthogonal  transformation;  therefore 

YE  =  ^ExAE  ='^$E^AE  =  -[®e^MEPAP 


Equating  the  two  above  expressions  for  Yg  results  in 


t  $e^mep  =  mep^®p^ 


or 


UE]  = 


mep^p^mpe 
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2.  3.  2.  3  Linearization  of  the  Differential  Equation 


The  linearisation  procedure  is  established  by  taking  the  partial  derivatives 

with  respect  to  each  error  source.  The  notation  used  is:  6  =  d/dt,  = 

th  ^ 

parti,  1  derivative  with  respect  to  the  i  error  source. 

Neglecting  the  higher-order  gravity  term,  we  find  that  the  navigation  system 
equation  to  be  linearized  is 

x=  --^5c  +  [MEpMpA]1[MApMpEJ2Xg  +  [M£a]  /  (A) 

XV 

When  the  following  assumptions  are  applied 

a.  p  is  a  known  constant 

b.  6Xg  =  0,  the  nominal  trajectory  reference 

c.  f  (A)  =  0,  products  of  small  perturbations  are  zero 

(see  Section  2 . 3 . 3 . 2  for  f  (A)) 

the  linearized  equation  becomes 

.  ~%65t  +  if  X6R  +  (M  a!M  +  6MEA1MAE}Xs  +  MEA6  /(X) 

R  R 

The  following  terms  in  the  above  expression  are  expanded  as 

/  2  2  2 
a.  R  =  J  X  +  Y  +  Z 

6r  =  XT6X/R 

b*  6MAE2  =  6^MAPMPE^2  =  6MAPMPE  +  MAP6MPE 


Since  ^MAp  is  in  error  due  to  a  misalignment  of  the 
accelerometer  true  input  axis  from  the  one  calibrated, 
it  is  treated  as  an  accelerometer  error  source.  Con¬ 
sequently,  6Map  =  0  for  these  equations.  The  matrix 
6Mp£  results  from  platform  angular  misalignments. 
Therefore 


6MAE2  =  MAP[6MPEJP  =  map{*pJmpe 


where 


[*p]  =  the  skew  symmetric  matrix  made  up  of  plat¬ 
form  angular  errors  expressed  in  platform 
coordinates 

C‘  6MEA1  =  6fMEPMPAJi  =  5MEPMPA  +  mep6mpa 

Since  Mp^  is  calculated  from  M^p,  there  is  no  error  in 
mPA.  that  is»6MpA  =  0.  Here  6Mep  results  from 
computer  errors.  Therefore 


6MEA1  =  tfiMEpV~PA 


=  -MEp[®p]M 


PA 


where 


[®p]  =  the  skew  symmetric  matrix  made  up  of 
computer  errors  expressed  in  platform 
coordinates. 


By  using  these  relationships,  the  equation  reduces  to 

rT 


6x  =  ”*3 


R  L  R 


3XX 
~Z - 1 


6X  +  iMEAMAp[*p]MpE 


-  MEp[®p]MpAMAE|xg  +  Mea  6 /(A) 
=  MC6X  +  MEp{®p  -®p}MpEX8  +  Mea6  /  (A) 


J  tj 


^(Symmetric) 


Using  the  properties  discussed  above,  these  equations  are  further  reduced  to 


6X  =  Mg6X  +  Ma  M£A6/  (A) 


where 


Xg,  Ygl  Zg  =  components  of  Xg 

4I  =  -  0  =  vector  expressed  in  ECI  coordinates,  which  is  the 

difference  between  computer  axes  and  platform 
axes,*  due  to  the  various  error  sources  (see 
Section  2.3.3) 

^  =  the  vector  expressed  in  ECI  coordinates  of  platform 
errors 

vi.  1 

This  notation  is  the  same  as  the  one  given  in  Reference  2,  o  161. 


0  =  the  vector  expressed  in  ECI  coordinates  of 
computer  errors 

6/"  (X)  =  the  vector  expressed  in  accelerometer  coordinates 
of  the  accelerometer  errors. 

These  then  constitute  the  basic  differential  equations  (variational  equations) 
for  computing  navigation  error  sensitivities,  as  a  function  of  the  navigation 
system  error  sources.  These  equations  are  put  into  a  pseudo-state  vector 
matrix  form  as  4 


x.  =  M  x.  +  F. 
i  xi  1 


where 


X1  X2  x3  are  position  partials  due  to  ith  error  source 

•  •  • 

x^  x,.  x^  are  XYZ  velocity  partials  due  to  i  error  source 

x7  x8  x9  are  orientation  partials  due  to  ith  error 

source  (also  see  Section  2.3.3) 
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o 


is  the  forcing  function  of  accelerometer  errors 
due  to  i*h  error  source 


forcing  function  of  orientation  rate  errors 
ith  error  source 


state  vector  representing  the  navigation 
sensitivity  due  to  i**1  error  source. 

2.3.3  Error  Sources 

The  error  sources  of  the  unaided  inertial  guidance  system  fall  into  the  general 
categories  of  initial  conditions,  accelerometei ,  gyro,  and  platform  errors . 

In  addition,  terminal  errors  are  included  to  be  able  to  account  for  guidance 
and  control  errors. 

The  general  notation  used  for  the  identification  of  error  sources  is  a  seven- 
character  symbol.  Not  all  are  explicitly  used  when  it  is  convenient  to  drop 
certain  characters  without  loss  of  generality,  or  when  an  error  source  is 
explicitly  denned  otherwise.  Each  character  is  defined  in  Table  1. 

Table  G-l  (in  Appendix  G)  gives  the  symbol,  description,  and  units  for  each 

$ 

error-source  type  presently  considered  .  Additions  to  this  list  are  easily 
accommodated  when  a  certain  component  does  not  fit  the  model  error  con¬ 
sidered  here.  For  purposes  of  describing  error-source  models,  the  phase 
index  is  dropped  and  the  error  source  is  represented  as  a  unit  vector  to 
derive  sensitivities,  that  is 


"1* 

"O' 

"O' 

EKfm;  EKI 1  = 

0 

,  EKI  ’  = 

1 

,  EKI  3  = 

n 

.0. 

_0_ 

1 

m  -J 

NOTE: 


One  of  the  I  indices  is  dropped  for  initial-  and  terminal -t,rpe  error 
sources,  since  I  »  3  for  these  categories. 


I 
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Table  1.  General  Notation  Used  for  Identification  of  Error  Sources 


Character 

1 

Symbol 

E 

Description 

Identifies  it  as  an  error  source  symbol 

2 

K 

Categorical  index  where  K  equals: 

I  -  initial  condition  error 

P  -  platform  error 

G  -  gyro  error 

A  -  accelerometer  error 

T  -  terminal  condition  error 

3  and  4 

1 

Numerical  ordering  of  the  error -source 
types  within  each  category 

1  =  00,  01,  02,  ....  99 

5 

m 

Identifies  a  component  or  axis  number 
m  =  1,  2,  or  3 

6  and  7 

n 

Phase  index  used  to  identify  at  what  time 
the  error  source  was  activated  (see  Sec¬ 
tion  2.3.4)  n  =  01,  ...  12 

(For  Example,  EA102-11  identifies  the  No. 2 
accelerometer  Type  10  error-source  active 
during  the  11th  phase  of  the  error  analysis.) 

2.3.3. 1  Initial  Condition  Errors 

As  indicated  in  Section  2.3.1,  there  are  two  options  fr  .  specifying  the  initial 
position  and  velocity  errors.  For  Option  1,  the  errors  cannot  be  explicitly 
defined  but  generally  would  be  the  same  as  for  Option  2  (with  a  change  in 
azimuth  direction).  In  the  definitions  for  initial -condition  errors  given  in 
Table  G-l,  it  is  assumed  that  Option  2  is  being  used  and  that  the  platform  is 
aligned  with  respect  to  vertical.  In  a  launch  from  an  earth  site  (tQ  £  0),  any 
uncertainty  of  the  launch  location  will  result  in  initial  -condition  errors  of 
position  and  velocity,  and  in  the  transformation  matrix  Mgp  of  the  computer. 
These  are  obtained  as  follows. 


2 . 3 . 3 . 1 . 1  Initial  Position  Error  s 


The  position  error  sensitivities  transform  into  ECI  coordinates  as 


(m  =  1,  2,  3) 


Initial  velocity  for  navigation  is  computed,  where  t  -  0,  as  follows 


V  =uXR  =  (w  1  )  x  (X  Y  +  Y_  Y  +  Z  1  ) 
o  o  e  u  c  u  °u  ou 


=  -«  Y  ”  1  w  X  Y  =  X  X  +  YY 

eou  eou  ou  ou 


where 


X  ,  Y  ,  I  =  ECI  unit  vectors 
u  u  u 

w  =  earth  rotation  rate 
e 

Therefore,  initial  position  errors  icoult  in  the  velocity  errors 


6X  =  -w  6Y  =  -w  x_ 
o  e  o  e  2 


=  x. 


6 Y  =  w  6X  - 
o  e  o 


«  X,  =  X 

el  5 
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\  ,*u 


or 


The  transformation  matrix  used  by  the  computer  is  developed  from  the  initial 
position  data,  therefore  errors  in  position  result  in  errors  in  Mj,p  {Mj-p  = 
Mp£  see  Section  2.3.1).  A  downrange  error  results  in  a  negative  rotation 
about  the  2-axis  of  and  a  cross-range  error  in  a  positive  rotation  about 

the  3 -axis.  There  is  no  error  due  to  altitude  errors.  The  angular  errors 
are  proportional  to  1/Rq;  therefore,  the  computer  error  in  launch  coordinates 
is 


This  is  transformed  into  the  ECI  coordinates 


When  t  >  0,  there  is  no  explicit  coupling  introduced  into  velocity  and 
o 

orientation  elements  of  the  state  vector  for  position  errors  and  the  elements 
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2.3.3. 1.2  Initial  Velocity  Errors 

Initial  velocity  errors  are  transformed  directly  into  ECI  coordinates  as 


/ 


6J 


*  ME1  [E!2m 


(m  =  1,  2.  3) 


X1X2X3X7X8X9  =  0 


2 . 3 . 3 . 1 . 3  Ip’tial  Platform -Orientation  Errors 

The  platform  errors  {^)  are  initial  rotations  about  the  platform  axes  and  are 
transformed  into  ECI  coordinates  as 


r 


-  41  -  ^Ep 


(m  =  1,  2,  3) 


x 


1  •  •  *  ‘ 


=  0 


2 . 3 . 3 . 2  accelerometer  Error  Sources 

In  the  design  and  manufacture  of  accelerometer  components,  every  attempt  is 
made  to  achieve  an  output  that  is  a  functio  .  only  of  input  acceleration  along 

one  of  its  axes.  Unfortunately,  this  is  never  achievable  due  to  inherent 

* 

d  sign  characteristics  and  manufacturing  tolerances.  The  general  equation 


See  Reference  2  for  a  discussion  of  design  characteristics,  etc. 


32- 


for  an  accelerometer's  output  at  a  given  time  can  be  written  as 

Ao  *  K0  +  K1A1  +  K2A?  +  K3A!  *  *4*2  +  Vs  +  K6A1A2  +  K7A1A3 
+  K8(A2  +  A3)1/2  +  Vl(A2  +  Al)1U  +  K10A2  +  KUA3 
+  K12A2A3  +  R 

where 

Aj  =  the  sensed  acceleration  along  the  defined  (theoretical) 
input  axis 

A^  =  sensed  accelerations  normal  to  A^ 

K j  =  coefficients  that  may  be  functions  of  time  and  are  the  result 
of  design  characteristics,  manufacturing  tolerances,  or 
environmental  effects  (vibration,  temperature,  a,...) 

R  =  a  remainder  term,  wmch  includes  all  higher -order  terms 
assumed  to  be  sufficiently  small  to  neglect. 

One  notable  remainder  term,  not  treated  here,  is  a  dynamic  response  to 

acceleration  transients.  Of  course,  any  one  component  design  does  not 

have  all  the  terms  presented  above;  and  only  those  that  are  significant  are 

included  in  a  given  error  analysis.  The  purpct_  oi  component  calibrations 

is  to  measure  these  coefficients  and  thereby  compensate  for  their  effects. 

♦ 

It  is  assumed  that  the  compensation  is  achieved  by  the  following  corrections 
to  the  accelerometer  output 


A  =  A 
oc  o 


(ka  +  k_  A2  +  . 

\  rtc  2c  o 


+  K_  A  -  + 
5c  o3 


Ki2cAo2Ac?) 


Compensations  could  also  be  achieved  by  biasing  the  target  conditions  in  the 
guidance  equations,  etc. 


where 


A  =  the  correcteu  output  of  an  accelerometer,  where  all  constants 

OC  and  accelerometer  outputs  have  been  scaled,  based  on  the 

calibrated  scale  factor  K,  of  each  accelerometer 

lc 

K  =  the  correction  coefficients  derived  from  instrument  calibration 
c  or  based  on  inherent  design  characteristics.  (Not  all  are 
determinable  from  the  above  methods.) 

A  .  =  (i  -  2,  3)  =  the  accelerations  normal  to  the  accelerome  c-r's 
01  input  axis,  derived  from  the  outputs  of  all  three  accelerometers 
as  follows 


o 


A 


o3 


=  M,  -A 

rs.  J  O 


where 

A  -  the  vector  of  accelerometer  outputs 
o 

A  _  =  a  vector  composed  of  the  acceleration  components  along  the 
0  2-axis  of  each  accelerometer 

A  ^  =  a  vector  composed  the  acceleration  components  along  the 
°  3 -axis  of  each  acc  ometer 

M  . and  M.0  are  defined  in  Section  2.3.1,  Since  the  A  coefficients  are 
n.7,  A3  1  c 

sufficiently  small  with  respect  to  Kj  (which  implies  linear  design),  their 
products  are  negligible  and  the  equation  reduces  to 


A 

oc 


(Kn-V'*KlAl  +  (K2-K2c)Al  + 


(K12  "  K12c)A2A3 
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1 


l  l 


Taking  the  partial  derivatives  of  this  equation  results  in 


6A  =  6Kn  +  5K.A.  +  6K,A,  +  .  .  .  .  6K.-A..A, 
oc  0  i  1  21  12  2  3 


where  6A.  =  0  (nominal  trajectory). 

Combining  the  equations  of  each  accelerometer  and  placing  them  into  the 
sensitivity  form,  the  final  vector  matrix  form  is 


AA11  ° 


6f (A)  =  6A  =  [  i]  (EAOOm)  +  0  AA21  0  {EAOlm} 


0  A, 


+  .  .  •  +  I  0 


AA12AA13 


AA22AA23 


{EM2m} 


m  =  (1,  2,  3) 


AA32AA33 


-/here 


A*,  —  A  .  „  ,  —  h'v  ,  r,X 

A1  A21  AE  s 


acceleration  components 
alon^  nput  axes 


/ 


AA22 


=  MA2XA1 


\ 


acceleration  components 
along  2 -axis 


A13 


A23 


LA 


A33  J 


=  MA3*A1 


acceleration  components 
along  3 -axis 


2. 3. 3. 3  Gyro  Erro„  Sources 

For  the  purpose  of  describing  the  error-source  model  for  gyro  components, 
it  is  convenient  to  discuss  the  model  in  terms  of  a  single -degree-of-freedom 
integrating  gyro.  The  general  equation  for  eyro  rates  can  be  written 

0  =  C  +  C.A  +  C_A  +  C , A  A  +  C.u,  +  C,u-  +  C,w,  +  C„A.,A, 

0  11  23  3  1  3  43  52  61  723 

+  C8A2  +  C9A1  +  C10A3  +  C11A1A2  +  R 


where 

Aj  A^  =  sensed  accelerations  along  the  input,  output,  and  spin 
reference  axes 

w.  w  ui  =  rates  about  the  input,  output,  and  spin  reference  ^.xes 

1  u  J 

=  coefficients,  which  may  be  functions  of  time  and  ai  ■  the 
result  of  design  characteristics,  manufacturing  tolerances, 
or  environmental  effects  (vibraucn,  temperature,  etc.) 

R  =  a  remainder  teem,  which  includes  all  higher-order  terms 
assumed  to  be  sufficiently  small  to  neglect. 
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One  notable  term  ia  dynamic  response  to  transient  inputs.  It  is  assumed  that 
the  platform  and/or  gyro  servo  loops  are  designed  with  sufficient  bandpass 

i 

and  static  gain  to  make  the  effects  of  transients  or  sustained  rotational 
dynamic  inputs  negligible. 

As  in  the  case  of  accelerometers,  not  all  coefficients  are  applicable  to  a 
given  design;  only  those  indicative  of  the  particular  components  are  considered 
in  any  one  analysis.  Also,  component  calibration  measures  some  of  these 
coefficients  and  the  compensation  for  their  effects  is  assumed  to  be  included 
in  the  navigation  system  equations.  The  compensation  method  assumed  is  the 

)Qt 

calculation  of  compensating  torquing  signals  io  the  gyros.  Thus,  the  com¬ 
pensated  gyro  rate  equation  is 


0  =  0  -  (Cn  +  C.  A  .  + 
c  Oc  lc  ol 


+  ^“oB  +  *  •  *  *  +  CilcAolAoE) 


where 


=  MGpMpAX0 

=  acceleration  along  gyro  input  axes 

Ao2 

=  MG2Xol 

=  acceleration  along  gyro  2-axes 

Ao3 

=  MG3Xol 

-  acceleration  along  gyro  3-axes 

Wol 

=  MGPWoP 

=  rates  about  gyro  input  axes 

Wo2 

=  MG2"ol 

=  rates  about  gyro  2 -axes 

Wo3 

=  MG3Wol 

=  rates  about  gyro  3 -axes 

A  is  as  defined  in  Section  2. 3. 3. 2 
o 


=  measured  or  computed  vector  rates  of  platform  axes 


M^..  is  as  defined  in  Section  2.3.1 


Use  of  the  compensation  of  tl  e  transformation  matrix  is  more  generally  in 
practice,  which  is  equivalent. 
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By  following  the  same  approach  a«  for  the  accelerometers,  the  final  equations 
for  gyro  error  aeneitivities  in  vector  matrix  form  become 


1  *  =  MEG?c 


where 


=  [I]  {EGOOm}  + 


Gil 

0 


G21 

0 


G3 1 J 


{EGOlm} 


u 


11 


0 

0 


CO. 


21 

0 


0 

0 


w. 


'31-» 


{EG  06  m) 


+  .  .  .  + 


AG11AG12 


0 

0 


AG  1AG22 


0 

0 


AG31AG32  J 


{EGllm} 


(m  =  1 ,  2,  3) 


where 


G1 


Gil 


G21 


G31 J 


= 

GE  8 


acceleration  components 
along  gyro  input  axes 
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acceleration  components 
along  gyro  2  -axis 


G12 


AG2L 


»-  "G32-I 


=  MG2SG1 


G13 


G23 


l-~G33 


=  MG3XG1 


acceleration  components 
along  gyro  3 -axis 


■»  m 

UJ 

11 

'  U1 ' 

W1 = 

U) 

21 

-MGP 

CO 

2 

rates  about  gyro  input  axes 

-31 J 

LW3- 

(a) 


12 


(0. 


22 


OJ 


»-  32-1 


=  MG2"l 


rate  components  about 
gyro  2 -axis 


00 


13 


go 


23 


00. 


33. 


r  V[ 

G3  1 


rate  components  about 
gyro  3  -axis 
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Li 

2 . 3 . 3 . 4  Platform  Errors 

In  addition  to  the  initial  condition  errors  discussed  in  Section  2.3.3. 1, 
acceleration-sensitive  errors  arise  due  to  the  structural  deformation  of  the 
gimbals  under  acceleration  loads,  and  to  static  servo  response  due  to  plat¬ 
form  mass  unbalances.  The  general  equation  for  platform  acceleration- 
sensitive  errors  can  be  written  in  the  vector  matrix  form 


*  ■ 

1 - 

o 

rsi 

Oe 

i 

> 

oo 

O 

_ _  _  1 

®2 

= 

AP3 

{EPOlm}  + 

Aw 

*3. 

.  o  Api. 

■  °  AP2 

fAP2AP3 

0 

1 

{EP02mJ 


L 


AP1AP3 


AP>AP2j 


f£P03m}  (m  =  1,  2,  3) 


V  = 


X 


9J 


=  MEp0p  (r\j,  -  0) 


where 


P 


P2 


= mpe*. 


L  AP3  J 


acceleration  't  -nponents 
along  platforrr  xes 
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2 . 3 . 3 . 5  Terminal  Condition  Er  ror  a 


The  terminal  condition  error*  transform  into  the  following  ECI  coordinate* 
by  utilizing  the  transformation  matrices  developed  in  Section  2.3.1. 

Position  Errors 


=  M£T1  lET.m)  (m  =  1,  2,  3) 


s  0 


Velocity  Errors 


MET2 {ET2m} 


(m  =  1,  2,  3) 


X1X2X3X7X8X9  =  ° 


2.3.4  Transition  Matrix 

Since  the  majority  of  error  sources  are  acceleration-sensitive,  their  forcing 
functions  (excluding  accelerometer  bias  and  gyro  bias  drift)  are  zero  during 
free -flight  sequences;  therefore,  it  is  more  efficient  computationally  to 
propagate  the  sensitivity  vectors  across  free  flight,  by  using  the  transition 
matrix  rathe ”  than  by  solving  each  error -source  sensitivity  independently. 
Additionally,  when  an  error -source  type  has  time -dependent  statistical 
characteristics,  it  becomes  both  convenient  and  efficient  to  subdivide  the 
total  trajectory  time  into  phases  and  treat  each  error  source  of  this  type  as 
an  independent  error  source  reinitialized  (x^(t)  =  0)  at  each  phase  time. 

Error  source  types  that  were  active  during  previous  phases  are  called 
inactive  vectors,  while  errors  sources  that  are  active  during  the  present 


phase  are  termed  active  vector*.  For  any  one  type  of  error  source,  all  it* 
inactive  vectors  are  updated  to  the  present  time  by  using  the  transition 
matrix(es)  and  combined  statistically  (see  Section  2.4.3)  to  derive  the  total 
effect  on  the  navigation  data  statistical  characteristics.  Thus,  the  following 
procedure  is  used  to  update  or  propagate  vector  sensitivities 

x.  (t)  =  ®(t,  t)x.(t) 


where  #  (t,  t)  is  the  transition  matrix  obtained  from  the  solution  of  the 
homogeneous  differential  equations 


#(t,  T)  =  M(tXf(t,  T), 


*(t,  t}  =  I 


where 


M(t)  = 


0  I  0 

M0(t)  0  MA(t) 

0  0  0 


M„  and  M.  are  defined  in  Section  2.3.2, 
G  A 


The  solution  of  this  equation  is  achieved  by  solving  the  homogeneous  differen 
tial  equations  (in  ECI  coordinates)  for  each  initial  condition.  Since  during 
free  flight  M^(t)  =  0,  the  solution  re  juires  only  six  independent  solutions; 
during  powered  flight  nine  independent  solutions  are  required.  In  addition, 
the  following  property  of  the  transition  matrix(es)  is  used 


«(t.t  )  =  ®(t.t.)«(t.t  ) 
i  k  i  j  j  k 
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To  make  an  error  analysis,  the  following  data  is  required  for  the  differential 
equations  and  the  error -source  equations: 


t 


time 


X  Y  Z  nominal  position  components 

in  ECI  coordinate  system 

•  •  • 

X  Y  Z  nominal  velocity  components 

in  ECI  coordinates 


X  Y  Z 

8  8  3 


nominal  sensed  acceleration 
components  in  ECI 


oj  |  u)  2  ^ 

MPE(tl 


nominal  platform  (body)  rates 
in  platform  (body)  coordinates 

direction  cosines  of  platform 
(body)  axes  with  respect  to 
ECI  axes 


This  data  is  generated  by  a  trajectory  program  for  a  particular  vehicle 
configuration  and  mission  requirement  and  is  written  on  a  magnetic  tape, 
which  constitutes  a  basic  input  for  the  error  analysis.  The  last  two  categories 
(u>p.  and  Mpg )  are  required  only  for  analyzing  a  strapped-down  configuration. 
For  a  torqued-platform  configuration,  the  rates  (“^  are  input  as  a  table  of 
rates  vs  time  and  Mp^,(t)  is  calculated  from 


M 


PE 


[Hp] 


M 


PE 


where  these  matrices  are  as  defined  in  Sections  2.3.1  and  2.3.2. 

When  error  propagation  must  be  evaluated  beyond  the  time  for  which  there 

is  data  from  the  trajectory  tape,  the  program  has  the  capability  to  generate 

the  required  data  [  M_  =  /(X  Y  Z  t) 3  for  calculating  the  transition  matrix  and 
^  G 


t 


drag  error  sensitivity.  The  data  is  generated  from  the  free-flight  equations 
of  motion,  based  on  a  simplified  oblate  earth  model  (see  Reference  5),  as 
follows 


where 


+  5*1(22:1  .  5*1  *  41 
Rr\?  7T  V 


“  z  J  A2  -  iL*l  z 


(.3 ♦ 12*)  --r'llt  «. 

\  hZ4-/  R2  NR2  7  /  1 


where  GM,J,A.  II,  and  D  are  nominal  constants,  defined  in  Table  G-2 
(in  Apoeudix  G). 

The  initia*  conditions  for  these  equation  are  ootained  from  the  trajectory 

* 

tape  x.  termination  e*  the  tape  data,  or  they  can  be  input.  Tne  three  criteria 


fui  terminating  the  equations  of  motion  are 
i  .  Specified  time  (t^,) 


An  option  is  available  to  terminate  (abort)  the  trajectory  tape  data  at  a 
turn*  prior  to  the  end  of  the  tape. 
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Specified  range  angle  from  termination  of 
the  trajectory  tape  data,  i.e. ,  when  6  = 

where 


8  =  coa 


XTX  +  YtY  +  ZTZ 


an<*  are  valuea  at  the  tape 
termination  time 

iNoie:  0  <  0T  <  v 

Specified  altitude  hT.  There  are  two  criteria  for  this 
termination:  The  trajectory  can  be  terminated  when 
h  -  h.p  and  the  slope  (h)  is  positive,  when  the  slope 
is  negative . 


DATA  PROCESSING  EQUATIONS 


2.  4 

2.  4.  1  Basic  Coordinate  Syttemi 

For  the  purpose  of  output  data  presentation,  two  coordinate  systems  are 
available:  the  ECi  coordinate  system  and  the  LH  (Local  Horizontal)  system. 
The  latter  is  sometimes  referred  to  as  the  orbit  plane,  and/or  the  Radial/ 
Tangential/Normai  (RTN)  coordinate  system.  The  transformation  between 
ECI  coordinates  and  local  coordinates  is  developed  from  the  nominal  trajec¬ 
tory  position  and  velocity  vectors.  In  local  coordinates,  X  is  defined  as 
down  ratige,  i.  e.  ,  as  directed  along  the  projection  of  the  inertial  velocity 
vector  onto  the  plane  normal  to  the  geocentric  radius  vector;  Y  is  vertical; 
and  Z  is  cross  range,  forming  a  right-hand  coordinate  system.  The  local 
coordinate  system  is  an  inertial  system  at  time  t.  defined  by  the  nominai 
conditions  only,  and  is  used  in  the  transformation  of  position,  velocity,  and 
orientation  sensitivity  vectors  as  well  as  those  of  transition  matrices  and 
covariance  matrices.  The  matrix  is  defined  as 


mel  =  mle 


wliere 


( 
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v  *  -VX  ♦  XV  t,  _  -Z(XX  4  y Y)  4  ZD2 
VE  C  ”n  "  - - RT5 - 


and  XYZXYZa  e  ECI  components  of  »hc  nominal  position  and  velocity 
vectors  at  time  t. 

NOTE:  4*  are  left-hand  rotations. 

It  is  also  convenient  to  form  the  matrix  as 

LE 


M 


LE 


“M 


LE 

0 


M 


LE 

0 


0 

0 


M 


LE  J 


mel  “  mle 


2.  4.  2  Vector  Errors 

Vector  errors  are  derived  by  scaling  the  sensitivity  vectors  by  an  appropriate 
constant.  The  constant  used  in  this  operation  is  given  the  symbol  <jv.  (i 
ranges  from  1  to  n,  which  is  the  total  number  of  error  sources,  i.  e.  ,  all 
active  and  inactive  vectors.  See  Section  2.  3.  4. )  The  units  of  a.  for  each 
type  of  error  source  are  given  in  Table  G-l.  Table  G-3  gives  the  conversion 
factors  used  by  the  program  for  scaling  <r .  into  the  units  of  the  sensitivity 
vectors.  The  operation,  therefore,  is 


AX.  =  r.KjXj 


where  the  last  expression  is  used  throughout,  implying  the  first. 


-47- 


AX.  ;•  the  error  vector  in  ECI  coordinates.  The  usual  implication  of  a  is 
that  it  represents  the  standard  deviation  of  the  particular  error  source. 
However,  it  could  also  represent  the  mean  value  of  an  error  source  or  be  a 
constant  that  produces  sensitivity  vector  output  in  any  desired  units.  When 
it  is  desired  to  output  the  vector  in  local  horizontal  coordinates,  the  following 
transformation  is  made 


LH 


2.  4.  3  Covariance  Matrix 

The  basic  equation  for  navigation  error  is 


"■£vi 

i=  1 

In  this  equation,  «.  is  the  magnitude  of  the  i^1  error  source.  Each  error 
source  is  considered  a  random  variable,  the  means,  variances,  <ir,d  correla¬ 
tions  of  which  are  assumed  to  be  known.  The  accuracy  of  any  navigation 
system  performance  is  measured  in  terms  of  the  probability  that  AX  (or  some 
function  of  AX  -  see  Section  2.  4.  5)  is  within  certain  specified  values.  To 
maximize  this  probability,  it  it  necessary  to  compensate  for  the  effects  of 
the  means  of  the  error  sources.  This  can  be  done  by  the  methods  discussed 
in  Section  2.  3.  3,  or  by  offsetting  the  guidance  constants  so  that  the  effects 
are  negated  at  some  specified  time  point,  for  the  latter  method,  the  <r  can 
be  input  to  represent  the  mean  value  of  each  error  source  and  the  resulting 
vectors  summed,  that  is 


n 

E(AX)  =  £-i4E(«i) 
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where 

E(«.>  =  <r-  representing 

*  4  * 

L  =  the  expectation  operator 

This  calculation  i*  not  made  in  the  program,  although  the  magnitude  of  the 
mean  (expected)  value  of  the  navigation  vector  can  be  obtained  from  the 
covariance  matrix  outpjit  described  below,  if  each  error  source  is  correlated 
with  the  other  by  unity.  Alternately,  the  error  vectors  can  be  summed  by 
utilizing  desk  calculators.  When  it  is  assumed  that  the  mean  values  have 
be&n  compensated,  the  second  statistical  moments  of  navigation  data  are 
determined  in  '-ovariance)  matrix  form  by 

=  E(x  ^  1  +  i2*2  +  ---)(x^€j  +  £*«  +  ---) 


t -it 

ECl  i=  1  j=  1  1  J  lJ  1  J 


where 


p..  =  the  correlation  coefficient  between  the  i  and  j  error 
J  source.  (In  this  way,  both  time  and  cross  correlation 
of  error  sources  are  handled. ) 

=  a  9  X  9  matrix  in  which  the  diagonal  elements  represent 
the  variances  of  the  navigation  data  and  the  off-diagonal 
elements  are  the  covariances. 

Since  this  is  a  symm  ,trical  matrix,  the  correlation  coefficients  of  the 
navigation  data  are  calculated  and  presented  in  the  lower  half  of  the  matrix. 
The  correlation  coefficient*  are  obtained  from 


mn 


<r  <r 
m  n 


m  ^  n 


-49- 


where  E(AX_AX  )  ie  the  mn^  elements  of  the  covariance  matrix  and  <r  , 
m  n  m 

•’  the  m  and  n  standard  deviations  of  the  m  and  n  coordinates  calcu¬ 

lated  from  the  square  too t  of  the  respective  diagonal  terms.  The  <r's  of  the 
navigation  data  are  also  presented  in  the  output. 

When  vector  errors  are  presented  in  LH  coordinates,  the  covariance  matrix 
is  also  presented  in  this  coordinate  system.  It  is  computed  by 


All  of  the  above  operations  presuppose  the  condition  of  linearity,  which  is 
usually  satisfied,  but  should  not  be  assumed  always  to  be  true.  (This  is 
mentioned  merely  as  a  precaution  about  the  underlying  assumptions  of  the 
program. ) 

Z.  4.  4  Transition  Matrix  and  Trajectory  Variables 

As  stated  in  Section  2.3  4,  the  transition  matrices  are  used  to  propagate  not. 
only  sensitivity  vectors  across  free  flight,  but  inactive  vectors  for  data  pro¬ 
cessing  as  well.  In  addition,  when  free-fligbt  time  histories  are  desired,  the 
transition  matrix  plays  a  paramount  role.  For  various  analytical  studies,  it 
is  desired  to  obtain  the  transition  matrix,  and  so  its  output  is  made  available. 
Like  the  previous  outputs,  it  is  available  in  either  the  ECI  or  LH  coordinate 
systems.  Since  it  is  computed  in  ECI  coordinates,  the  operation  required  to 
present  it  in  local  coordinates  is  easily  obtained  from 

AX.(t)  =  *(t,  t)aK.(t) 


AXt(t) 


LH 


MLE(t)A*i(t)  = 


=  Mle  (t)*(t.  t)Mel(t)43?.(t) 


®LH*t,T)  =  Mle^*^  t)mEL(t) 


LH 
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In  addition  to  the  transition  matrix,  it  is  convenient  and  sometimes  necessary 
to  know  the  reference  trajectory  conditions  for  which  the  transition  matrix  is 
valid.  These  are  presented  at  both  times  (t  and  t)  and  in  both  ECI  and  local 
reference  coordinates,  the  latter  being  defined  and  computed  from  ECI 
coordinates  as 


LAT  (deg) 


Geocentric  latitude 


sin 


-l  Z 


2<lat<2 


LONG  (deg)  Longitude  measured  positively  east  from  Greenwich 
=  tan  *  +  0^  •’  «gt  0  £  LONG  <  2tt 

where  0^  is  used  to  reference  the  ECI  coordinates  to  Greenwich  when  the 
trajectory  reference  is  not. 

ALT  (n  mi)  Altitude  in.  nautical  miles  above  the  surface  of  an 
oblate  earth 

P  -  Re 

=  £0767  TW5 


where 


R  = 

e 


_ A(1  -  e) 

(1  +  (e2  -  2e)  cosZ\)1/2 


See  Section  2.  3.  5  for  definitions. 


VEL  (ft/ sec)  Inertial  velocity  magnitude 


•  2  -2  -2  1/2 
=  (X  +  Y  *  Z  )  ' 


FPA  (deg)  Flight  path  angle,  defined  as  the  angle  the  inertial 
velocity  vector  makes  with  the  local  geocentric 
horizontal 


._-l  XX  +  YY  +  ZZ 
=  sin  - - gv - 


-  \  <  FPA  <  J 
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AZ  (deg)  Azimuth  of  inertial  velocity  vector  measured  clockwise 
from  north. 


=  tan 


-1  VE 


0  <  AZ  <  2ir 


See  Section  2.4.  1  for  definition  of  V_  and  Vkt. 

E  N 

2.  4.  5  Mission  Evaluation 

This  is  an  optional  output  of  the  program  and  operates  on  the  final  or  end 
conditions  of  a  particular  error  analysis.  It  is  called  upon  when  a  specified 
mission  termination  criterion  is  to  be  evaluated.  The  parameters  necessary 
to  evaluate  mission  success  can  occur  in  a  variety  of  categories,  only  a  few 
of  which  are  considered  here. 

2.  4.  5.  1  Fixed  Altitude 

This  criterion  is  used  primarily  for  re-entry  evaluation  and  presents  the 
coordinates  of  down- range  (MD)  and  cross- range  (MC)  miss  at  a  fixed  altitude 
with  respect  to  earth  fixed-target  coordinates;  in  addition,  the  time  dispersion 
(MT)  is  presented.  The  orientation  of  the  down- range  miss  coordinate  is 
along  the  projection  of  the  relative  velocity  vector  onto  the  plane  normal  to 
the  target  geocentric  radius  vector.  The  method  of  computing  these 
quantities  is 


VR  =  V  -  VT 

=  (X  +  w  Y)X  +(Y-uX)Y  +ZZ 
e  u  e  u  u 

=  Xj  +  Yn  Y  +  Z„Z 
R  u  R  u  R  u 


where: 


X  Y  Z 


u  u  u 


XYZXYZ 


X 


R 


Y 


R 


=  target  velocity  in  ECI  coordinates 
=  relative  velocity  in  ECI  coon,  nates 
=  ECI  unit  vectors 

=  components  oi  the  nominal  position  and  velocity 
vectors  c.t  the  target 

=  tne  components  of  VD  in  ECI  coordinates 
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The  relative  velocity  vector  i«  transformed  into  the  coordinate  system  at  the 
target  by 


y  «i 

*RR 

*xr" 

*RR 

=  MRE1 

*r 

L  ^'RR* 

*^R  - 

where 


M 


REi 


0 

O-b 

S4* 


0 

-S4j 

r.< 


C\ 

9 

i*-SX 


JL 


o 

i 

0 


sx 

0 

cx 


~C0  S0 
-50  C0 
0  0 


S0  =  ^  SX  =  £  = 


z 


V 


RE 


RH 


C0  =  w  CX  =■£  C*  =  'RIN 


V 


RH 


D  =  7?+  Y* 


YXr  +  XYr 


RE 


RN 


-Z(XRX  +  YYr)  +  ZrD 
- - - - 


VRH  =  >/VRE  "  VRN 


and  X  Y  Z  XR  YR  ZR  are  ECI  components  of  the  nominal  target  coordinates 
and  relative  velocity  vector  at  time  t. 
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Thun,  the  relative  velocity  components  in  the  target  coordinates  are:  XDD 

KR 

the  altitude  rate,  YDD  the  velocity  component  in  the  deiined  cross -range 
direction  (i  zero),  and  Zpp  in  the  defined  down- range  direction.  When 
V0_  =  =  0,  vertical  re-entry,  4*  Is  undefined.  In  tnis  case,  ip  Is  set 

to  ze?o  so  that  down-range  displacements  would  be  directed  north  and  cross¬ 
range  east. 


Based  or.  the  above  definitions  for  a  coordinate  system,  a  position  error 
vectoircan  be  transformed  into  this  coordinate  system  to  determine  the 
altitude,  cross-range  and  down-range  dispersions  at  the  nominal  time  by 


■AXt‘ 

"AX- 

AYt 

mrei 

AY 

t4 

N 

<1 

_ 1 

AZ 

where  AXj  AY,j,  AZ^.  are  altitude,  cross-range,  and  down- range  dispersions 
at  the  nominal  time  and  AX  AY  AZ  are  ECI  dispersions  for  a  given  error 
source. 

To  derive  the  first-order  dispersions  for  the  condition  of  zero  altitude  error 
(AX^,  =  U),  the  following  constraint  equation  is  csed 

AXT(t  +  At)  =  0  =  AXT(t)  +  XRR(t)At 


from  which  the  time  dispersion  is  calculated  as 


T 


axt 

XF„R 


At  =  M 


T 


Using  this  relation,  the  cross-range  and  down-range  misses  at  time  t  +  At 
are  obtained  by 


AYT(t  +  At)  =  Mc  =  AYT(t)  +  tRR(t)At 


^  AY 


T 


AZT(t  +  At)  =  Md  =  AZT(t)  +  ZRR(t)At 

..  A„  zrr  aY 
-  AZ_  t  ax~, 

XRR 

Thus 


r  * 

“ 

*M1 

0 

n 

-axt- 

» 

X 

to 

to 

Mc 

= 

0 

1 

0 

ayt 

LmdJ 

ZRR 
.  *RR 

0 

1 

- 

_AZ  x_ 

"mt" 

'AX' 

Mc 

=  M„. 

■L 

R1MRE1 

AY 

.md. 

_AZ_ 

The  position  vector  error  for  each  error  source  is  thus  transformed  into  its 
dispersion  parameters  at  the  target.  The  covariance  matrix  of  these 
parameters  is  developed  oy  partitioning  the  ECI  covariance  matrix  to  include 
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I 


only  the  position  variances  and  covariances  from  which  the  following 
covariance  matrix  is  determined 


£■ 


P 


=  MTR1MRE1  2Hmtrimrei) 


where 


t  - 

ME1 


the  covariance  matrix  for  mission  evaluation  Option  1 


t- 


-  the  partitioned  covariance  matrix  in  ECI  coordinates 


ECI 

P 


The  nomir  1  trajectory  conditions  are  also  listed  with  this  output  and  include 

LAT,  LONG,  ALT  As  described  in  Section  Z.  2.4 
VEL/R  (ft/ sec)  Magnitude  of  relative  velocity  ector 


FPA/R  (deg) 


*■«  r  +  *r  +  4 


Relative  flight  path  angle,  defined  as  the  angle  the 
relative  velocity  vector  makes  with  the  local  horizontal 


=  3in 


XXR  +  YYr  +  ZZ 


TO 


—  -  £  <  FPA/R  <y 


AZ/R  (deg) 


Azimuth  of  the  relative  velocity  vector,  measured 
clockwise  from  north 


V 

=  tan'1  0  <  AZ/R 

VRN 


/ 

V 


2.  4.  5.  2  Fixed-  range  Angle 

This  criterion  is  similar  to  the  fixed-altitude  case,  except  the  relative  range 
coordinate  (AZ -)  is  fixed  and  the  altitude,  cross-range,  and  time  dispersions 
are  determined.  The  constraint  equation  is 

AZT(t  +  At)  =  0  =  AZx(t)  +  ZRR(t)At 


from  which  the  time,  cross- range,  and  altitude  dispersions  are 


•■•t  =  M 


T  " 


AYT(t  +  M)  =  Mc  =  AYt 
AXT(t  +  At)  =  My  =  AXT(t)  +  XRR(t)At 

XRR 

=4XT-r*“T 


Thus 


r*  •m 

“ 

1 

m 

m  m 

Mt 

0 

n 

axt 

AX 

yj 

z 

Mc 

= 

0 

i 

0 

AYt 

"  MTR2MRE1 

AY 

M 

i 

o 

*RR 

N 

AZ 

•  m 

V 

^RR 

.. 
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.... 

Vov-tf/vvA-,  -v- >./v3- 


The  vector  error  output#  and  covariance  are  handled  in  the  same  way  as  in 
the  fixed"  altitude  case. 

2.  4.  5.  3  Generalised  Linear  Transformation 

«ntwTMnMM*Ni — w— M*wuii  i— W  jiuub 

This  option  is  used  when  linear  transformations  between  ECI  position  and 
velocity  errors,  and  some  other  parameters  are  known,  e.g.  ,  the  midcourse 
maneuver  velocity  sensitivities,  as  a  function  of  injection  errors  for  a  space 
probe.  A. matrix  is  developed  from  (up  to  10)  input  matrices,  and  then  used 
to  transform  injection  errors  into  the  desired  parameter  errors.  Each  error 
source  vector  is  transformed  and  a  covariance  matrix  is  calculated  from  the 
partitioned  ECI  covariance  matrix,  computed  as  follows. 

The  matrix  is  formed  from 


M  =  [Mj  [M^]  •  •  •  [Mj  n<  10 

The  vector  errors  are  transform*  i  by 


m  m 

1 

’ax' 

2 

AY 

3 

*  [M] 

uZ 

4 

AX 

5 

AY 

6 

J 

A Z 

The  covariance  matrix  is  calculated  as 


2t  =m£mT 

ME3  ECI 
PV 
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2.  4.  6  Platform  Reference  Attitude 

When  running  a  torqued  platform  or  strapped-down  case,  reference  platform 
orientation  is  automatically  presented  as  a  function  of  time.  Listed  in  the 
output  are  the  first  two  rows  of  the  platform  matrix  (Mp^,)  and  three  angles 
defined  and  computed  as 

THETA  (0  in  deg)  Pitch:  the  angle  that  the  platform  l-axis  makes 

with  the  local  horizontal  plane 

+  l <m  +  zaz,)..sesw 

where  (IX),  ( 1 Y) ,  (1Z)  are  the  elements  of  the 
first  row  of  Mp^  representing  the  direction 
cosines  of  Pj  with  respect  to  the  ECI 
coordinates  axes 

PSI  (4>  in  deg)  Heading  or  Yaw:  the  angle  that  the  projection 

of  the  platform  1-axis  onto  the  horizontal 
plane  makes  with  north,  positive  clockwise 

=  tan*  1  0  £  4/  <  2ir 

where  lc 

,*r  _  -z[X(lX)  +  Y(  1 Y)]  +  DZ(1Z) 

IN - EE - 

PHI  (0  in  deg)  Roll:  the  angle  that  the  platform  2-axis  makes 

with  the  local  horizontal  plane 

where  (2X),  (2Y),  (2Z)  are  the  elements  of  the 
second  row  of  M  ^  representing  the  direction 
cosines  of  with  respect  to  the  ECI 
coordinates  axes. 
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SECTION  3 


COMPUTER  PROGRAM  INPUT/OUTPUT 

3.  1  INTRODUCTION 

This  section  is  intended  as  a  guide  to  users  of  the  error  analysis  program 
in  preparing  the  input  data.  Brief  descriptions  of  the  program  capabilities 
and  the  procedures  for  data  input  are  presented,  and  the  available  output 
data  formats  are  described.  There  are  standard  forms  for  preparing  the 
necessary  data  input,  axxd  they  are  shown  in  Appendix  A. 

The  program  actually  consists  of  two  separate  programs:  The  first,  called 
ERAN,  computes  error  sensitivities  based  on  given  trajectory  data,  naviga¬ 
tion  system  instrument  configuration,  and  a  model  of  the  component  error 
sources.  Trajectory  data  is  supplied  via  a  tape  generated  by  a  trajectory 
program(s)  (II- STAGE,  TRIP,  MVS,  etc.).  The  instrument  configuration  is 
supplied  as  the  input  data  of  component  orientations.  The  error  models  are 
supplied  on  an  error- source  schedule  sheet  that  specifies  the  error  «ources 
to  be  considered  for  this  configuration  and  when  they  are  active.  Based  on 
this  information,  the  ERAN  program  computes  sensitivity  coefficients,  and 
outputs  these  onto  an  ERAN  tape. 

The  second  program,  called  0UTP,  basically  processes  the  ERAN  tape  into 
the  prescribed  formats.  Inputs  to  this  program  are  the  logical  controls, 
i.  e. ,  the  desired  frequency  of  output,  coordinate  systems,  formats,  etc.  , 
and  the  specification  of  standard  deviations  of  the  error  sources,  along 
with  any  correlation  coefficients  between  error  sources.  The  ERAN  tape 
can  be  processed  many  times  without  rerunning  the  ERAN  Program. 


3.2  ERAN  INPUT  DATA 

3.  2. 1  Trajectory  Tape 

Prior  to  running  the  error  analysis  program,  a  trajectory  tape  must  be 
prepared  that  contains  the  following  information  relative  to  the  nominal 
mission: 

Time 

Nominal  position  coordinates  in  the  ECI  system 

Magnitude  of  position  (radius)  vector 

Nominal  velocity  coordinates  in  the  ECI  system 

Nominal  sensed  acceleration  coordinates  in  the 
ECI  system 

Rates  of  platform  (body)  axis  in  platform  (body) 
coordinates  (optional) 

Direction  cosines  of  the  platform  (body) 

axes  with  respect  tc  the  ECI 

coordinate  system  (optional) 

Each  file  cn  the  trajectory  tape  contains  the  information  for  one  ERAN  case 
and  consists  of  the  following  logical  records: 

1 st  Record 

Word  1  ~  1B27 
Word  2  =  N 

Word  3  =  0  (N  =  siz  ;  of  a  data  record  i.  e.  ,  11  or  23) 


t 

x,  y,  z, 

R 

X  Y  Z 


W1  w3 

IX  1Y  1Z 
2X  2Y  2Z 
3X  3Y  3Z 


Continuous  i>ats  Point  Record 


Word  1  -  1 1335 
Word  2  =  t 
Word  3  =  X 

Word  4  =  Y 


Word  5  =  Z 
Word  6  =  R 
Word  7  =  X 

s 

Word  8  =  Y 

s 
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Word  9 

=  Z 

8 

Word  17  =  1Y 

Word  10 

=  X 

Word  18  =  1Z 

Word  11 

=  Y 

Word  19  =  2X 

Word  12 

=  Z 

Word  20  =  2Y 

Word  13 

*  W1 

Word  21  =  2Z 

Word  14 

sw2 

Word  22  =  3X 

Word  15 

=  w3 

Word  23  =  3Y 

Word  16 

=  IX 

Word  24  =  3Z 

(Words  13  to  24  are  omitted  if  N  =  11) 

Left  Side  of  a  Trajectory  Discontinuity 
Word  1  =  3B35 

Words  2  through  24  are  the  same  for  all 
data  records 

Right  Side  of  a  Trajectory  Discontinuity 
Word  1  =  5B35 

Words  2  through  24  are  the  same  for  all 
data  records 

Last  Record  of  a  Trajectory 

Word  1  =  1B29 
EOF 

All  Aerospace  Corporation  trajectory  programs  are  mechanized  to  prepare 
a  trajectory  tape  in  the  proper  format  for  input  to  ERAN.  The  tape  writing 
intervals  should  be  no  greater  tnan  4  seconds  during  powered- flight  phases 
and  32  seconds  during  coast  or  free-flight  phases.  Higher  tape  densi¬ 
ties  would  cause  no  problem,  but  lower  densities  would  tend  to  degrade  the 
accuracy  of  the  integrations.  The  integration  step  size  used  for  ERAN  is 
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not  necessarily  the  game  as  that  of  the  input  tape  interval,  but  is  controlled 
by  input  data.  The  ERAN  Program  interpolates  between  trajectory  data 
points  to  obtain  proper  values  for  integration.  The  integration  routine  used 
by  ERAN  is  based  on  a  fourth- order  Runge-Kutta  method. 

3.  2.  2  Error  Sources 

The  Error  Source  Schedule  (see  Table  A-l  in  Appendix  A)  is  used  to  identify 
which  of  the  available  error  sources  are  to  be  run,  and  for  which  time 
periods  (phases)  they  are  to  be  considered.  Table  G- 1  defines  the  symbol 
and  units  for  each  error  source  in  the  order  it  appears  on  the  Error  Source 
Schedule.  Initial  and  terminal  error  sources  are  listed  individually,  while 
component  (accelerometers  and  gyros)  and  platform  errors  are  listed  as 
errc  types.  Thus,  when  a  component  or  platform- error  type  is  considered, 
sensitivities  for  all  three  components  or  axes  are  automatically  and  indepen¬ 
dently  run.  Each  source  or  type  of  error  can  be  considered  in  one  or  all  of 
12  possible  independent  phases  of  the  trajectory.  This  phase  capability  is 
provided  to  accommodate  time- correlated  errors  and/or  independent  error- 
source  magnitude  changes.#  It  becomes  increasingly  mre  important  as 
the  mission  time  duration  increases. 

Sections  3.  2.  2.  1  through  3.  2.  2.  5  present  a  summary  of  the  error- source 
types  and/or  error-  model  equations  (see  Section  2.  3.  3). 

3.  2.  2.  1  Initial  Condition  Errors 

El  11  -*  El  13  Initial  Position  Errors 

EI21  -*EI23  Initial  Velocity  Errors 

EI31  -*  EI33  Initial  Platform  Alignment  Errors 

If  the  initial  time  of  the  trajectory  is  equal  to  or  less  than  zero,  the  pro¬ 
gram  assumes  a  launch  from  an  earth  fixed  pad  and  automatically  calculates 


Phase  logic  will  be  used  to  identify  and  control  "Reset,  "  a  feature  that 
updates  navigation  data  as  a  function  of  external  measurements  (see 
Appendix  D). 
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an  initial  velocity  and  vertical  alignment  error  consistent  with  the  initial 
position  error. 

3.  2.  2.  2  Accelerometer  Errors 


AA  =  EAOO  +  EAOKAj)  +  EAOZ^2)  +  EAOS^3) 

+  EA04(A2)  +  EA05(A3)  +  EA06(A1A2) 

+  EA07(AjA3)  +  EA08(AZ2  +  A32)1/2 
+  EAC9(Aj)  (A,2  +  A32)1/2  +  EA10(A22) 
+  EAll(A32)  +  EA12(A2A3) 


whare 

AA  =  the  i1*1  accelerometer  error  (i  =  1,  2,  3) 

EAi  =  (f  =  00,  01,  ....),  as  defined  in  Table  G-l 

th 

Aj  =  input  axis  acceleration  of  i  accelerometer 
A^  =  acceleration  components  normal  to  input  axis 

3.  2.  2.  3  Gyro  Errors 

0  =  EG00  +  EG01(A1)  +  EG02(A3)  +  EG03(A1A3) 

+  EG04(w3)  +  EG05(u>2)  +  EG06(uj) 

+  EG07{A2A3)  +  EG08(A2)  +  EGO^Aj2) 

+  EG1C^A32)  f  egh(a1a2) 
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where 


til 

0  =  the  i  gyro  drift  error  (i  =  1,  2,  3) 


3.2.  2.4 


EG f  =  (f  =  00,  01, 


.),  as  defined  in  Table  0-1 


th  tli 

A.  =  acceleration  along  j  axis  of  i  gyro 
tli  tli 

w.  =•  rat.e  about  j  axis  of  i  gyro 

J 


j  =  1  input  axis 
j  =  2  output  axis 
j  =  3  spin  reference  axis 
Platform  Errors 


typical 

identification 


0.  =  EPOli(Aj)  +  EP02i(Ak)  +  EP03i(A^Ak) 


where 


.th 


0.  =  platform  angular  error  about  i  axis  (i  =  1,  2,  3) 

th 

Aj,  Ak  =  platform  acceleration  components  norma?  to  t  axi« 
EPf  =  (f  =  00,  01,  ....),  as  defined  in  Table  G- 1 


3.  2.  2.  5  Terminal  Errors 

ETll-°-ET13  Terminal  Position  Errors 

ET21  — •  ET23  Terminal  Velocity  Errors 


Note:  These  errors  can  be  considered 

only  on  the  last  phase  of  the  error 
analysis  run  and  can  be  applied 
only  at  the  trajectory  tape  abort 
time  (see  Section  3.2.  3)  or  at  the 
end  of  the  trajectory  ti  pe. 


The  procedure  for  filling  out  the  error- source  schedule  is  as  follows: 

a.  Determine  which  error  sources  are  to  be  considered. 


fa.  Establish  if  arty  of  these  error  sources  are  to  have  non¬ 
unity  autocorrelation  functions;  i.  e. ,  establish  whether 
multiphase  logic  is  required, 

c»  In  the  first  column,  insert  an  "X'1  in  each  row  element 
that  describes  the  error  soui.  e  to  be  considered, 

d.  In  the  second  (and  following)  columns,  insert  an  "X"  in 
the  row  elements  corresponding  to  the  error  sources 
that  require  phase  logic  control  for  that  particular  phase* 
(see  example). 

The  error  sources  that  have  been  called  for  by  inserting  X's  in  the  appro¬ 
priate  squares  in  Column  1  will  be  initialized  and  become  active  at  the  start 
of  the  case,  TSUBO  (see  Section  3.2.  3.  5).  Insertions  of  X's  into  the 
appropriate  squares  in  Column  2  will  cause  those  error  sources  to  be  re¬ 
initialized  at  the  start  of  the  second  phase  (i.  e.  ,  when  the  time  of  the 
simulation  is  equal  to  the  value  entered  into  TCOP  (see  Section  3.2.  3.  5). 

The  sensitivity  vectors  from  the  first  phase,  for  those  error  sources  will 
then  become  inactive.  These  inactive  vectors  will  be  updated  at  the  desired 
output  times  by  the  transition  matrix  and  be  combined  statistically  with  the 
active  vectors  to  derive  the  total  effect  on  the  navigation  data  statistical 
characteristics  (see  Section  2.  3.4).  Similarly,  X's  in  the  third  column 
will  cause  those  error  sources  to  be  reinitialized  when  the  time  of  the 
simulation  is  equal  to  the  value  entered  into  "^GOP  +  1,  etc.  In  this  manner, 
up  to  12  independent  sensitivity  vectors  can  be  ere  ted  for  each  error 
source  on  the  schedule  sheet 

This  completes  the  input  on  the  Error  Source  Schedule  sheet.  It  essentially 
indicates  to  ERAN  which  error  source  sensitivities  tu  calculate  and  whether 
any  of  these  types  require  phase  logic  control. 

o 

It  will  be  necessary  in  setting  up  the  input  data  for  QKJTP  *0  assign  sigma 
values  (the  standard  deviations  of  each  error  source)  and  correlation 

*Error  sources  that  change  standard  deviations  or  have  time-varying 
correlation  coefficients. 


coefficient#,  when  applicable,  between  error  sources.  The  identification 
of  a  given  sigma  value  with  a  sensitivity  vector  is  accomplished  by  mentally 
assigning  a  number  to  each  error  source.  The  numerical  ordering  of  error 
sources  is  determined  by  starting  in  the  first  column  of  the  Error  Source 
Schedule  and  counting  down,  then  going  to  the  second  column,  and  so  forth. 
Note  tliat  three  independent  error  sources  are  associated  with  each  "X"  in 
an  element  of  a  component  or  platform  error, 

3.  2.  3  Orientation  and  Control  Data 

The  data  sheet  used  to  set  up  the  platform  and  component  orientations  and 
to  input  the  necessary  program  control  data  and  constants  ;s  shown  as 
Table  A-2  in  Appendix  A.  Control  data  and  constants  have  preassigned 
values;  therefore,  only  those  numbers  that  deviate  from  them  need  be  input. 
The  symbols  used  are  summarized  in  Table  G-2,  along  with  their  pre¬ 
assigned  numerical  values  and  units. 


3.  2.  3.  1  Initial  Platform  Alignment 

There  are  two  options  for  specifying  the  initial  platform  orientation  (see 

Figure  2).  The  first  assumes  that  the  platform  1-axis  (F^)  is  aligned  along 

the  geocentric  vertical.  The  azimuth  orientation  is  specified  through  an 

input  of  PSIP  (4 >  ,  a  left-hand  rotation  in  platform  coordinates).  With  nc 
P 

input  in  PSIP,  the  platform  would  be  aligned  so  that  the  3-axis  (P^)  would 

be  north  and  the  2-axis  (P~)  east.  With  a  positive  ip  ,  the  P-.  axis  would  be 

c  pi 

rotated  toward  the  east 


The  second  option  allows  for  a  more  general  initial  platform  orientation; 
here,  three  angles  are  specified  for  aligning  the  platform.  The  initial 
alignment  is  such  that  the  P, ,  P^,  and  axes  are  along  the  ECI  X,  Y,  and 
Z  axes,  respectively.  PHIP  (0  )  rotates  the  platform  positively  about  its 
3-axis;  LAMP  (X.^)  rotates  the  platform  negatively  about  its  2-axis;  and 
P,*>IP  (iji  )  rotates  it  negatively  about  its  1-axis,  in  that  order.  The  option 
is  used  when  it  is  desired  to  align  to  geodetic  cr  astronomic  latitude  as  a 
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vertical  reference  for  ground  alignment,  or  when  platforms  are  assumed 
to  be  aligned  in  orbit  with  some  stellar  instruments. 

3.  2.  3. 2  Initial  Conditions 

There  are  two  options  for  initial  condition  specifications  (see  Figure  5). 

If  no  input  to  PSII  (4>t)  given,  the  program  assumes  the  initial  conditions 
to  be  referenced  to  platform  axes;  thus,  down  range  is  along  the  axis, 
cross  range  along  the  P^  axis,  and  altitude  along  the  axis.  If  an  input 
of  PSII  (^)  is  specified,  the  program  assumes  that  vertical  or  altitude 
errors  are  along  the  geocentric  vertical,  and  down- range  errors  are 
referenced  from  north.  Cross- range  errors  are  therefore  (fcj  +  90*  from 
north. 

3.  2.  3.  3  Gyro  Orientation 

The  initial  orientation  of  the  gyros  and  their  axes  is  illustrated  in  Figure 
Gyro  alignment  is  made  by  specification  of  an  axis  of  rotation  (1,  2,  or  3) 
and  an  argument  (angle)  of  rotation.  The  program  allows  up  to  five 
independent  rotations  in  any  order  desired.  Each  rotation  operates  on  the 
gyros  as  a  triad;  i.  e. ,  all  three  gyros  are  being  rotated  and  thus  maintain 
their  axis  orientation  with  respect  to  each  other  during  these  rotations. 

The  axes  of  rotation  referred  to  above  are  those  of  the  No.  1  gyro.  Upon 
completion  of  this  set  of  rotations,  there  remains  an  additional  degree  of 
rotational  freedom  of  each  gyro  about  its  input  axis,  specified  by  PSI.  (ty.) 

(i  =  No.  1,  2,  or  3  gyro). 

3.2.  3.  4  Accelerometer  Orientation 

There  are  two  options  available  for  the  alignment  of  accelerometers  (see 
Figure  4).  The  first  is  the  specification  of  an  orthogorai  triad  and  the 
method  is  identical  with  that  described  for  the  gyro  components;  i.  e. ,  an 
axis  and  argument  are  specified  for  aligning  the  accelerometer  input  axes. 
Then  an  additional  degree  of  freedom  about  each  accelerometer's  input  axis 
is  specified  by  BETA^  ((3^)  (i  =  No.  1,  2,  or  3  accelerometer). 
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The  second  option  allows  for  nonorthogonal  accelerometer  configurations. 

Mere  the  method  of  specification  (of  an  axis  and  an  argument)  is  the  same; 
however,  each  accelerometer  is  specified  independently.  The  initial 
orientation  of  each  accelerometer  is  the  same  with  its  1  - ,  and  3-axes 
along  those  of  platforms  Pj,  an<*  ^3*  respectively.  The  data  sheet  for 
nonorthogonal  accelerometers  is  shown  as  Table  A- 3  in  Appendix  A.  To 
exercise  this  option,  a  non-zero  entry  must  be  made  in  data  location  field 
ACCEL  10;  conversely,  a  zero  entry  negates  the  option. 

3,  2.  3.  5  ERAN  Control  Data 

(JUT  This  entry  controls  the  ERAN  tape  writing  frequency.  If  no 

entry  is  made,  the  tape  writing  frequency  for  this  case  will  be 
two  records  per  trajectory  discontinuity  and  two  records  per 
ERAN  phase  discontinuity.  If  a  non- zero  entry  is  made,  a 
single  record  will  be  written  at  aach  multiple  of  100  sec  for 
powered  flight,  and  of  1000  sec  for  free  flight,  unless  other¬ 
wise  specified  in  PPF  and  PFF. 

PPF  This  entry  controls  the  tape  writing  interval  during  powered 

flight  for  values  other  than  the  standard  100  sec,  when  0UT  *  0. 

PFF  Similarly  this  entry  controls  the  tape  writing  intervals  during 

free  flight  for  values  other  than  the  standard  1000  sec. 

TSUBO  This  is  initial  time  and  can  be  any  time  greater  than  or  equal  to 
the  first  time  point  on  the  trajectory  tape  (file).  If  the  first  time 
point  (t)  on  the  tape  is  not  zero  (t  =  0),  then  the  desired  starting 
time  must  be  entered. 

TSUBA  This  is  abort  time  and  can  be  any  time  less  than  the  last  .ime 
point  on  the  tape  (file).  If  an  entry  is  omitted,  the  tape  will  be 
processed  from  TSUBO  to  the  end  of  the  trajectory  tape  (fil^). 

TRAJ  When  the  trajectory  tape  contains  more  than  one  trajectory,  this 

entry  identifies  the  trajectory  (file)  to  be  processed  (e.  g.  ,  if  the 
entry  is  N,  ERAN  will  process  the  N*h  file  on  the  tape).  If  omit¬ 
ted,  the  next  trajectory  will  be  processed.  For  the  first  ERAN 
case  this  would  be  File  1. 

When  consecutive  files  are  being  processed  on  a  trajectory  tape, 
starting  with  File  1,  the  program  will  run  most  efficiently  when 
no  entry  is  made  to  TRAJ.  When  an  N  entry  is  made,  the  program 
will  process  the  N  file  for  that  case  and  all  subsequent  cases  until 
TRAJ  is  altered  by  input. 

ENDC  This  entry  controls  the  use  of  the  equations  of  motio  .  in  ERAN. 

These  equations  model  an  oblate  atmosphere -free  earth,  used 


1 

MAXT 

BTNP 

DTNF 

BMT 

BRTAB 


to  propagate  errors  beyond  abort  time  TSUBA.  The  pair  of 
entries  (ENBC  and  1)  control  the  termination  of  the  propagation. 

The  options  for  inputting  to  ENDC  are  as  follows; 

terminate  at  abort  time 

terminate  at  that  value  of  time  (sec),  which 
is  specified  in  the  next  entry  (>TSUBA) 

terminate  at  the  value  of  range  angle  (deg) 
beyond  the  termination  of  the  trajectory  tape, 
which  is  specified  in  the  next  entry 
(0  <  9  <  180) 

terminate  at  that  value  of  altitude  (ft)  with  a 
position  slope,  which  is  specified  in  the  next 
entry 

same  as  (d)  above,  with  negative  slope. 

This  identifies  the  location  field  where  the  numerical  value  of  the 
termination  control  is  to  be  input. 

This  entry  controls  the  maximum  running  time  of  the  equations  of 
motion.  It  is  preset  to  36,  000  sec;  i.  e. ,  when  t  =  36,  000  the  run 
will  be  aborted  unless  a  greater  value  is  entered. 

This  entry  controls  the  ER  AN  integration  step  size  during  powered 
flight.  No  entry  will  cause  the  program  to  integrate  at  its  nominal 
4-sec  integration  step.  * 

This  entry  controls  the  ERAN  integration  step  size  during  free 
flight.  No  entry  will  cause  the  program  to  integrate  at  its 
nominal  32 -sec  integration  step.  * 

This  is  the  flag  used  to  tdentify  a  case  where  the  platform(body)axes 
will  not  be  inertially  oriented  during  the  run.  A  non- zero  entry 
will  cause  the  program  to  seek  one  of  the  options  described  below. 

This  flag  is  used  to  identify  the  option  to  be  used  for  obtaining  plat¬ 
form  turning  rates  and  platform  direction  cosines.  A  non-zero 
entry  will  cause  the  program  to  determine  platform  rates  from  a 
tahle  of  input  rates.  From  this  data  the  direction  cosines  are 
derived  by  integrating  the  matrix  differential  equation  of  direction 
cosines.  No  entry  will  cause  the  program  to  read  this  data  (rates 
and  direction  cosines)  from  the  trajectory  tape. 


a.  No  Entry 

b.  TIME 

c.  THETA 

d.  ALTP 

e.  ALTM 


$ 

Note:  The  program  integration  routine  converges  on  each  multiple  of  the 
tape  writing  interval  when  0UT  0.  Therefore,  when  the  value  of 
PPF  is  less  than  DTNP,  the  former  would  be  the  integration  step 
size  used  in  powered  flight.  Similarly,  when  the  value  of  PFF  is 
less  than  DTNF,  the  former  is  used  for  the  free  flight  integration 
step  size. 
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TQOP 


3. 2. 3. 6 

0MEG2 
A 

GM 
e 
J 
H 
D 

MU 

The  numerical  values  of  these  constants  are  given  in  Table  G-2. 

3. 2.  3.  7  ERAN  Case  Control  Data 

Since  multiple  casen  from  one  trajectory  tape  can  be  run  sequentially  in 
using  the  ERAN  program,  two  cards  are  necessary  to  instruct  the  program 
as  follows: 

END  The  preceeding  cards  contain  all  the  data  necessary  to  run 

this  case. 

ENDJ0B  This  is  the  last  case  processed  by  ERAN.  Since  it  is  preprinted 
on  the  standard  form,  it  must  be  crossed  out  for  all  cases 
except  the  last. 

3. 2. 4  Tabular  Input 

3.  2.4.  1  Turning-  rate  Table 

The  standard  form  for  platform  turning  rates,  which  the  program  uses  if 
the  3RTAB  flag  (see  Section  3.  2.  3)  is  non- zero,  is  shown  as  Table  A-4  in 
Appendix  A.  The  definitions  of  symbols  and  the  method  to  be  used  to  input 
data  are  as  follows: 


This  entry  and  the  ten  that  follcw  it  are  used  to  identify  the 
time  to  terminate  a  phase.  No  entry  is  needed  to  terminate 
the  last  phases  consequently,  for  cases  in  which  there  is  only 
one  phase,  no  entry  is  required. 

\ 

Earth  Model  Constants 

— f—  — — —  ■  ■ '  — — ■ 

rotation  rate  of  the  earth 

equatorial  radius  of  the  earth 

gravity  constant  used  in  the  equations  of  motion 

ellipticity  of  the  earth 

constant  in  the  earth's  potential  function 

constant  in  the  earth's  potential  function 

constant  in  the  earth's  potential  function 

gravity  constant  (equals  GM)  used  in  the  variational  equations. 


ORDER 


refers  to  the  order  of  data  interpolation  to  be  used  by  the  pro¬ 
gram  to  establish  rates  between  data  inputs.  A  .1  entry  will 


cause  the  program  to  use  linear  interpolation,  a  &  quadratic, 
etc.  The  interpolation  routine  used  is  a  kth  order  Lagrangian 

identifies  the  total  number  of  time  points  in  the  table  that 
follows 

first  time  point  of  table  (t^  £  TSUBO) 


last  (N  j  time  point  of  table  (t^  £  TSUBA) 
rate  abour  platform  (body)  1  -axis  at  time  tj 


rate  about  platform  (body)  1-axis  at  time 
rate  about  platform(body)2-axis  at  time  tj 


rate  about  platform  (body)  2 -axis  at  time 
rate  about  platform(body)3-axis  at  time  t^ 


rate  about  platform(body)  3-axis  at  time  t^ 

Note:  If  the  table  contains  many  zeros  in  sequence,  they  can 
be  entered  by  writing  a  Z  in  the  prefix  field  and  the 
number  of  zeros  to  be  generated  in  the  value  field.  As 
an  example,  the  table  for  a  case  of  20  time  points,  zero 
rates  about  the  1-  and  3-axis,  and  the  first  two  rates 
about  the  2- axis,  also  zero,  would  look  like  Table  2. 
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Table  2.  Sample  for  a  Case  of  20  Time  Points 


PRJC  DOC 


Valve 


Remarks 


value  of  t. 


value  of  t. 


vaxne  of  w2(t3) 


20  zero  rates  about  1-axis 
and  2  zero  rates  about 
2- axis 


value  A 
20 


20  zero  rates  about  3- axis 


Also  note  that  the  reverse  side  of  the  standard  form  can  be  used  to  continue 
the  rate  table. 

3.  2.  4.  2  Equation  of  Motion  Initialization 

The  program  is  mechanized  so  that  the  equations  of  motion  can  be  initialized 
independent  of  a  trajectory  tape  input.  This  feature  is  used  when  it  is 
desired  to  obtain  transition  matrices  or  to  use  one  of  the  mission  evaluation 
options  to  derive  miss  coefficients.  The  format*  for  this  data  is  given  as 
Table  A-  5  in  Appendix  A. 


A  printed  standard  form  is  not  available. 


✓ 


3.  2.  j  Multiple  Cases 

As  mentioned  previously.  ERAN  has  the  capability  to  run  multiple  cases 
from  a  single  trajectory  tape.  The  data  used  for  the  first  run  is  retained 
for  the  second  (and  subsequent)  runs;  thus,  only  data  that  requires  changes 
from  the  preceding  runs  needs  be  input.  When  it  is  desired  to  eliminate 
the  effects  of  an  orientation  option  used  in  a  previous  case,  it  is  necessary 
to  input  a  minus  zerc  in  an  appropriate  location.  The  three  options  end  the 
methods  of  cancelation  are  as  follows: 


a. 

Platform  Orientation 
Option  2 

input  minus  0  in  PHIP 

b. 

Initial  Condition 
Orientation  Option  2 

input  minus  0  in  PSH 

c. 

Nonortho  go  nal 

input  minus  0  in  ACCEL  10, 

Accelerometer 

which  is  the  input  location  for 
the  first  axis  of  rotation  of  the 

No.  2  accelerometer  component 

This  last  operation  will  negate  the  logic  that  was  set  up  by  the  previous  non- 
orthogonal  case  and  will  therefore  interpret  the  data  for  the  No.  1  com¬ 
ponent  as  that  required  for  a  triad. 

Extreme  caution  should  be  exercised  when  attempting  to  change  the  Error 
Source  Schedule  for  an  operation  where  there  are  more  than  f>  phases. 

Some  knowledge  of  the  input  routine  is  necessary  to  present  the  intrinsic 
problem. 

The  D  option  (i.  e.  ,  D  in  the  prefix  field  of  the  inpu*  word),  used  to  input  the 
Error  Source  Schedule,  causes  two  words  to  be  stored  in  the  computer,  with 
the  last  6  characters  being  stored  in  the  location  immediately  following  the 
location  of  the  first  6.  When  there  are  no  entries  in  Columns  7  through  12, 
however,  only  the  first  6  characters  are  stored  and  the  second  location 
remains  unaltered.  It  is  then  apparent  that  an  X,  entered  beyond  Column  6 
for  a  previous  case,  cannot  be  eliminated  without  entering  at  least  one  X 
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Mm*  other  column  beyond  the  6th  for  the  case  in  question.  Cancella- 
ol the  error  source  for  all  phases  can  be  achic.ed  by  entering  zeros 
the  appropriate  locations  for  the  sigma  value  (see  Sections  3.  3.  3. 1  and 
Phase  logic,  for  an  error  source  during  phases  7  through  11,  can  be 
eliminated  when  there  are  less  than  12  phases  to  the  case  by  entering  an  X 
in  Cohrac  12.  Should  that  be  the  only  entry  on  the  line,  the  error  source 
would  he  eliminated  for  the  entire  case. 


OUTPUT  DATA 


3.  3 

The  output  data  processor  program  (0UT P)  takes  the  data  generated  by  ERAN 
and  produces  output  data  at  the  required  times,  with  the  prescribed  trans¬ 
formations  and  the  proper  format.  The  options  available  for  the  above  data 
follow. 

3.  3.  1  Output  (Print)  Times 

Option  0  Output  the  data  only  at  the  terminal  condition  of  the  case. 

Option  1  Output  at  the  phase  discontinuities  plus  the  terminal  conditions. 

Option  2  Output  data  called  for  by  Options  0  and  1  and  at  all  trajectory 

discontinuities  where  the  sensed  acceleration  goes  from  non-zero 
to  zero  or  from  zero  to  non-zero,  including  the  initial  time. 

Option  3  Process  every  time  point  on  ERAN  tape  as  determined  by  the  tape 
density  control. 


3.3.2  Output  Coordinate  Systems 


ECI 


LH 


EVALU 


Presents  data  in  an  Earth  Centered  Inertial  System,  where  the 
Z-axis  is  the  earth's  polar  axis  and  the  X-  and  Y-axes  are  in  the 
equatorial  plane.  Generally,  the  convention  is  that  the  X-axis 
passes  through  the  Greenwich  meridian  at  time  zero,  and  Y  com¬ 
pletes  a  right-hand  system;  however,  these  coordinates  are 
determined  by  the  particular  trajectory  program  used  to  generate 
the  input  tape. 

Presents  data  in  a  local  horizontal  coordinate  system,  which  is 
inertial  and  developed  from  the  nominal  trajectory  position  and 
velocity  vectors.  X  is  down  range,  i.  e.  ,  directed  along  the 
projection  of  the  inertial  velocity  vector  onto  the  plane  normal  to 
the  radius  vector;  Y  is  vertical,  i.  e.  ,  along  the  geocentric  radiu6 
vector;  and  Z  is  cross  range,  forming  a  right-hand  coordinate 
system. 

Presents  additional  data  at  the  terminal  condition  only  with  a 
prescribed  transformation.  *  Presently  there  are  the  following 
three  options  for  this  output: 

EVALU1:  Presents  the  down-range  (Mp),  cross-range  (M^), 

and  timing  (M-p)  errors  (misses)  at  a  fixed  altitude 
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A  special  format  is  used  for  these  transformations  (see  Section  3.  3.  3.  4). 


Presents  the  cross-range  (M^),  vertical  (My),  and 
timing  (M<p)  errors  (misses)  at  a  fixed  range 

General  —  represents  transformation  developed  from 
input  matrices. 


3.3.3 


sut  Data  Formats 


The  five  present  formats  for  output  data  will  be  described  in  Sections  3.  3.  3.  1 
through  3.  3.  3.  5.  Examples  of  the  formats  are  presented  with  the  test  cases 


i 

I 


in  Appendix  C» 


3.  3.  3.  1  Vector  Errors 


In  order  to  present  the  vector  errors  at  a  particular  time,  0UTP  first  update 
all  inactive  vector  sensitivities  (previous  phase  sensitivity  vectors)  by  pre- 
muitiplying  them  by  an  appropriate  transition  matrix.  Next  it  scales  the 
vector  sensitivities  by  the  proper  sigma  level  (input  data  of  the  sigma  of  a 
particular  error  source).  This  results  in  a  vector  error  in  the  ECI  coordinal 
system.  Finally,  it  performs  a  coordinate  transformation,  when  required. 
This  data  is  presented  in  a  standard  format  where 


Line  1  gives  the  run  date  and  job  identification  (see  Section  3. 4) 

Line  2  is  the  case  identification 

Line  3  is  the  nominal  time  and  coordinate  Bystem  identification 

Line  4  is  column  headings  where  DPX,  DPY,  DPZ  are  delta-position 

coordinates  to  th*»  nearest  ft 

DVX,  DVY,  DVZ  are  delta-velocity  coordinates  to  the  nearest 
0.01  ft/sec 

DOX,  DOY,  DOZ  are  delta -platform  orientation  coordinates  to 
the  nearest  0.  1  Bee 

Line  5-up  is  vector  error,  where  the  left  column  identifies  the  vector  as 
fQ  ilow  s: 


xxxx 


XX 


Four  characters  are  used 
to  identify  the  error -source 
type,  and  the  component  or 
axis  it  represents.  The 
identification  generally 
follows  the  symbols  given 
in  Table  G-l  with  the  fol¬ 
lowing  changes: 

Initial  Condition 

Accelerometers 

and  gyros 

Platform 

Terminal  Conditions 


Two  characters  are  used 
to  identify  the  phase  in 
which  this  particular  er¬ 
ror  source  was  initiated 
(01  to  12) 


£  is  replaced  by  O 

E  in  dropped  and  the  error  type  is 
followed  by  1,  2,  or  3,  indicating 
which  component  it  represents 

E  is  dropped  and  the  error  type  is 
followed  by  1,  2,  or  3  indicating 
which  platform  axis  it  represents 

E  is  replaced  by  T,  and  T  is 
replaced  by  O 


Note  that  the  order  in  which  the  error  vectors  are  presented  is  that  given  in 
Section  3,  2.  2. 


3  3,3.2  Covariance  Matrix 

The  covariance  matrix  is  formed  from  the  expression 


n  n 


ECI  i=  1  j=  1 


p.  .cr.«r.x.x. 
1  J  1  J 


T 


where 


X 

ECI 


=  covariance  matrix  in  ECI  coordinates  (9X9  matrix) 


n 


=  total  number  of  error  sources 


pij 


=  correlation  coefficient  between  the  itft  and  j**1  error  source. 
When  i=j,  p=l;  when  i^j,  p=0  unless  input  otherwise 


<r.:o%  =  standard  dwviationa  of  the  i^1  and  j^‘  error  sources 

»  sensitivity  vector  of  x*1  error  source 

— T  th 

x  j  =  transpose  of  j  sensitivity  vector 


To  obtain  the  covariance  matrix  ir,  the  local  horizontal  coordinate  system,  the 
following  operation  is  performed 


=  covariance  matrix  in  local  horizontal  coordinates 

=  matrix  (9  X  9)  which  transforms  a  sensitivity  vector 
from  ECI  to  local  coordinates 

=  transpose  of 

The  presentation  of  this  data  at  a  pa*ucular  time  is  in  a  standard  format, 
where 

Line  1  =  case  identification 

Lins  C  -  nominal  time  and  coordinate  system  identification 

Line  3  identities  it  as  the  covariance  matrix 

Line  4  =  column  heading  (same  definition  as  vector  errors) 

.Lines  5  to 

13  =  covariance  matrix  in  the  following  format: 

diagonal  presents  variances  of  navigation  errors  in  floating  point 

upper  elements  present  covariances  of  navigation  errors  in  float¬ 
ing  point 

lower  elements  present  correlation  coefficients  in  fixed  point. 


where 


% 

LH 


M 


LE 


M 


LE 


ftsmo 


I 


Line  14 


Line  15 
line  16 


3.  3.  3.  3 


=  standard  deviations  (sigmas)  of  the  navigation  errors  formed 
from  the  square  root  of  the  variances 

=  trajectory  variables  in  EC1  coordinates 

=  trajectory  variables  in  earth  reference  coordinate  system, 
where: 

LAT  =  geocentric  latitude  (deg) 

LONG  =  longitude  from  Greenwich  meridian  (deg) 

ALT  =  altitude  (n  mi)  above  the  surface  of  an  oblate 
earth  (n  mi) 

VEL  =  inertial  .  locity  (ft/sec) 

FPA  =  flight  path  angle  defined  as  the  angle  the  inertial 
velocity  vector  makes  with  the  local  geocentric 
horizontal  (deg) 

AZ  =  azimuth  of  the  inertial  velocity  vector,  measured 
clockwise  from  north. 

Transition  Matrix 


The  transition  matrix  is  used  in  the  ERAN  Program  to  propagate  sensitivity 
vectors  across  the  free-flight  sections  of  a  trajectory,  rather  than  integrating 
each  set  of  error-source  equations.  *  It  is  used  by  0UTP  when  presenting  a 
time  history  of  vector  errors  during  free  flight  in  the  same  manner,  i.  e.  ,  to 
propagate  the  sensitivity  vectors.  0UTP  also  uses  the  transition  matrix  to 
update  an  inactive  vector  (one  generated  in  a  previous  phase),  when  running 
multiphase  cases. 

The  transition  matrix  is  generated  in  the  usual  manner  -  b”  solving  the 
homogeneous  differential  equations  (in  ECI  coordinates)  for  each  initial  condi¬ 
tion  error.  To  obt«*>  the  transition  matrix  in  the  local  horizontal  coordinate 
system,  the  following  operation  is  performed: 


*(t,r)  =  M.  (tWt.rjM.  —  l(r) 
LH  —  ECI 


jg 

In  the  present  formulation  of  the  program,  it  is  assumed  that  the  accelerom¬ 
eters  are  disconnected  at  termination  of  a  powered  phase. 


wh«r« 


G 


i  -< 


•(t,  v) 

LH 


=  transition  matrix  in  local  horizontal  coordinates 


#(C,  t)  *  transition  matrix  in  ECI  coordinates 

ECI 


=  matrix  at  time  t,  which  transforms  a  sensitivity 
vector  from  ECI  to  local  coordinates 


M 


jjj,  T(t)  =  transpose  of  MT  Tr  at  time 


There  is  an  option  to  control  the  output  of  the  transition  matrix  (see  Section 
3.4).  When  called  for,  the  matrix  will  be  presented  at  each  discontinuity 
(phase  or  trajectory  tape).  In  addition,  if  a  time  history  (print  Option  3) 
is  called  for,  the  transition  matrix  will  be  presented  at  the  same  times  as 
the  vector  errors  and  covariance  matrix.  The  presentation  of  the  transition 
matrix  is  in  a  standard  format,  where 


Line  1 
Line  2 


Line  3 
Line  4 
Line  5 
Line  6-14 
Line  15 
Line  16 
Line  17 


Line  18 


3.  3.  3.4 


=  case  identification 

=  identifies  it  as  the  transition  matrix  and  the  applicable  time 
arguments  (t,  t) 

=  the  coordinate  system  identification 
=  format  identification 
=  column  headings 
=  transition  matrix 

=  trajectory  variables  in  ECI  coordinates  at  time  t 

=  trajectory  variables  in  ECI  coordinates  at  time  t 

=  trajectory  variables  in  earth  reference  coordinate  system 
at  time  t 

=  trajectory  variables  in  earth  reference  coordinate  system  at 
time  t 

Mi 8 8 ion  Evaluation 


As  indicated  in  Section  3.  3.  2,  there  are  presently  three  options  available  for 
presenting  data  at  the  end  of  the  case  that  can  be  used  for  evaluating  the 

$  T  - 

Except  for  free  flight,  a  transition  matrix  is  not  computed  in  Phase  1. 


z'1 

o 
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success  of  the  mission.  The  date  presented  when  one  of  these  options  is 
called  for  consists  of  vector  errors  and  a  covariance  matrix.  The  format 
for  each  option  is  as  follows: 


EVALU1 


Presents  data  at  the  reference  altitude  instead  of  the  reference 
time 


Vector  Errors 


Line  1 
Line  2 
Line  3 


Line  4 -up 


case  identification 

nominal  time  and  criterion  (ALT) 

column  headings,  where: 

MT  »  timing  error  to  the  nearest  0,  001  sec 

MC  =  cross-range  miss  tc  the  nearest  ft 

MO  =  down- range  miss  to  the  nearest  ft 

vector  errors,  where  the  left  column  identifies 
(described  in  Section  3.  3.  i)  the  vectors 


Covariance  Matrix 


Line  1 
Line  Z 
Line  3-5 


Line  6 
Line  7 


Line  8 


identifies  it  as  a  covariance  matrix 
column  heading 

covariance  matrix  output  (same  format  as  described 
in  Section  3.  3.  2) 

sigma  values 

headings  for  nominal  trajectory  conditions  at 
termination 

trajectory  conditions  where: 

LAT,  LONG,  ALT  are  as  defined  in  Section  3.  3.  3.  Z 

VEL/P  =  magnitude  of  relative*  velocity  vector 
(ft/sec) 

FFA/R  =  magnitude  of  relative  flight  path  angle  (deg) 
AZ/R  =  &*imuth  of  relative  velocity  vector  (deg) 


Velocity  vector  with  respect  to  rotating  earth. 


BVALU2  Presents  data  at  the  reference  range  angle  instead  of  at  the 
-  reference  time.  The  data  format  is  the  same  as  for  EVAJLUl 
except  as  noted  below. 

Line  2  criterion  (ALT  replaced  by  RANGE) 

Line  3  MD  is  replaced  by  MV  (vertical  miss  to  the  nearest  ft) 

EVALU3  Thle  option  is  used  for  special  cases  in  which  the  linear  transfer* 
(nation  between  ECJ  position  and  velocity  errors  and  some  arbi- 
trary  parameters  are  known;  e.  g. ,  the  midcourse  maneuver 
velocity  components  as  a  function  of  injection  errors  for  a  space 
probe,  some  orbit  elements,  etc. 

Vector  Error* 


Line  1 
Line  2 

Line  3 
line  4-up 


case  identification 

nominal  time  and  identification  of  EVALUATION 
OPTION  3 

column  heading  (1,  Z  ...  .  6) 

vector  errors  (floating  point)  where  the  left  column 
identifies  (as  described  in  Section  3.  3.  3.  1)  the  vector 


Covariance  Matrix 
-  -  —  —  ■■■■■«■■» 


3.  3.  i.Z. 


Line  1 
Line  2 
Line  3-8 

Line  9 


covariance  matrix  identification 
column  headings  (1,  Z  ...  .  6) 

covariance  matrix  (same  format  as  described  in 
Section  3.  3.  3.  2) 

si<gma  values 


Platform  Reference  Attitude  Time  History 


When  a  torqued  olatform  or  strapped-down  case  is  being  run,  a  time  history 
of  attitudes  and  direction  co  ines  are  given  at  the  end  of  the  case.  The  times 
are  the  same  as  those  called  for  by  the  print  option  The  format  for  presenta¬ 
tion  of  this  data  is  as  follows: 


Line  1  identification  of  type  of  output 
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Line  2  column  headings,  where 


THETA  (0) 
PSI  (i|i) 

PHI  <*) 


=  angle  (deg)  platform  1-axis  makes  with  the  local 
horizontal  plane 

=  angle  (deg)  the  projection  that  the  platform  )-axis 
onto  the  horizontal  plane  makes  with  north 

-  angle  (deg)  the  platform  2 -axis  makes  with  the 
local  horizontal  plane 


For  a  strapped-down  case  these  angles  would  be 
missile  pitch,  yaw,  and  roll  angles,  respectively. 

IX,  1Y,  12  =  direction  cosines  of  the  platform  1-axis,  with 
respect  to  the  ECI  coordinate  system 

2X,  2'Y,  22  =  direction  cosines  of  the  platform  2-axis,  with 
respect  to  the  ECI  coordinate  system 

Line  3 -up  time  history  of  above  data 


L 
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3.4  0UTP  INPUT  DATA 

Tbm  itudkrd  form  for  input  data  to  0UTP  ia  presented  as  Table  A-6  in 


Appendix  A. 
follows: 


Th«  definitions  and  procedure  for  filling  out  this  sheet  are  as 


HJ0BI  and  IU0B2 
FCA8E1  and  HCA5E2 


CASE 


S0PT 

PHIL 

EVALU 

F0PT 

SIGMA 


Up  to  40  characters  (30  each),  which  will  be  printed 
as  the  first  line  of  vector  error  output 

Up  to  60  characters  (30  each),  which  will  be  printed 
as  the  second  line  of  vector  errors  and  first  line  of 
all  other  outputs 

Print  time  option  as  discussed  in  Section  3.  3.  1.  No 
input  will  produce  Option  O. 

Controls  output  coordinate  system.  No  input  results 
in  LH  coordinate  system.  +ECI  results  in  both  ECI 
and  LH  output.  -ECI  outputs  only  in  ECI  coordinates- 

The  number  of  the  case  (N)  to  be  processed,  N  being 
the  N**1  trajectory  processed  by  ERAN  (but  not 
necessarily  the  Nth  trajectory  on  trajectory  tape). 

When  outputting  consecutive  ERAN  cases,  starting 
with  Case  1,  the  program  will  operate  most  effectively 
if  no  entry  is  made  to  CASE.  When  an  N  entry  is 
made,  the  program  will  process  the  N^  case  until 
CASE  is  altered  by  input. 

A  non-zero  entry  will  result  in  storing  the  sigma 
values  for  the  next  (and  subsequent)  runs;  whereas  a 
zero  entry  will  result  in  setting  all  sigma  entries  to 
zero  after  the  present  case  has  been  completed. 

Correction  term  for  longitude  output.  To  be  used 
when  the  trajectory  ECI  system  is  not  referenced  to 
Greenwich. 

Identifies  the  mission  evaluation  option  (if  any)  to  be 
output. 

Noi.\-*ero  entry  will  result  in  transition  matrlx(es) 
output. 

In  the  LOC.  field,  the  vector  number  is  input  and  the 
sigma  value  for  that  vector  is  put  in  the  value  field.  * 
Only  when  a  SIGMA  value  changer  is  it  required  to 
input  a  new  value;  otherwise,  it  will  assume  the  sigma 


^  '  """ 

An  entry  of  zero  will  cause  that  vector  to  be  eliminated  from  the  case. 


i 


RH0 


END 

ENDJ0B 


value  of  the  previous  vector  error.  As  an  example, 
if  unit  aensitivites  are  desired,  then  a  1  in  the  IOC. 
field  and  a  1  in  the  VALUE  field  will  result  in  output 
with  unit  scaling  of  all  error  sources.  *  However, 
when  the  S0PT  option  is  used  for  multiple  cases,  all 
desired  changes  to  the  sigma  table  must  be  entered 
explicity.  In  the  above  example,  all  desired  changes 
from  their  assigned  unit  values  wouin  have  to  be 
entered;  e.  g. ,  a  change  of  the  sigma  value  for  the 
first  error  source  would  only  alter  that  value,  all 
others  retaining  their  unit  values. 

These  are  correlation  coefficients.  Two  entries  are 
required  to  input  a  correlation  coefficient:  The  first 
identifies  the  error  sources  (by  vector  number)  that 
are  correlated,  and  the  next  gives  the  value  of  the 
correlation  coefficients.  All  correlation  coefficient 
data  is  retained  in  storage;  therefore,  when  running 
multiple  cases,  care  must  be  exercised  to  not  get 
unwanted  correlation  into  the  covariance  matrix 
calculations.  A  double-aero  in  the  field  for  assigning 
vector  numbers  of  a  correlation  coefficient  will  result 
in  eliminating  the  effects  (if  any)  of  that  previously 
stored  correlation  coefficient,  as  well  as  of  all  those 
that  followed  on  the  input  sheet. 

Same  control  as  discussed  f  r  ERAN 

Same  control  as  discussed  for  ERAN 


0UTP  has  the  same  logic  of  data  storage  (except  as  noted  in  the  S0PT  option) 
for  multiple  cases  as  ERAN.  Thsretcre,  only  chi. ages  to  data  need  be 
entered  for  runs  following  the  first.  To  negate  the  SYSTM  optJon,  six  seros 
must  be  entered  (i.  e. ,  when  ?  H  output  alone  is  desired  after  some  other 
option  on  the  previous  case  has  been  chosen). 


Note  units  of  error  sources  in  Table  C-i. 
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When  the  EVALU3  option  i«  used,  the  format  for  the  input  data  sheet  is 
protested  as  Table  A -7  ia  Appendix  A.  The  linear  transformation  matrix  [  M] 
used  i«  this  option  is  formed  from  products  of  input  matrices  by 

[M]  =  [MjHIAjHMj]  •••  [Mn] 
n  <10 

where  is  a  (6  X  6)  u.atrix  input. 


A  preprint  standard  form  is  not  available. 

( 
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SECTION  4 

SAMPLE  CASES 


Three  teat  cases  were  designed  to  demonstrate  the  procedures  for  filling 
out  input  data  sheets  when  exercising  the  various  program  options  available, 
and  to  present  data  in  all  of  the  output  formats.  The  data  sheets  used  to  set 
up  the  test  cases  are  presented  in  Appendix  B  and  the  output  listings  from 
these  runs  in  Appendix  C. 


The  trajectory  used  for  the  test  cases  was  one  that  had  been  designed  ko  place 
a  payload  into  a  24-hour  synchronous  equatorial  orbit.  Following  are  the 
major  trajectory  sequences: 


0  -  464 

464  -  477 
477  -  498 
498  -  1380 
1380  -  1685 
1685  -  20177 

20177  -  20288 


Powered  flight  from  launch  to  booster 
burnout 

Separation  sequence  (coast) 

Inject  into  100-mile  parking  orbit 

Coast  to  first  equatorial  crossing 

Inject  into  transfer  orbit 

Coast  in  transfer  orbit  to  apogee  of 
19, 300  n  mi 

Inject  into  synchronous  equatorial  orbit 


The  salient  features  of  the  test  cases  and  the  pertinent  input/output  options 
used  to  obtain  these  features  are  now  described. 


4.1 


TEST  CASE  1 


TW*  fa*  tfh*  «vak*hoi>  of  the  uncertainty  of  instantaneous  impact  prediction 
CUP)  for  ymmitttre  thrust  termination,  fay  using  the  inertial  navigator  data 
and  aswniag  a  vacuum  re-entry.  The  same  option  could  be  used  for  eval¬ 
uating  faefHsrtic  missile  accuracy  or  guided  re-entry  vehicle  accuracy  at  a 
fixed  altitude. 

The  configuration  of  the  inertial  navigator  was  one  in  which  the  platform 
1-axis  was  aligned  with  the  geodetic  vertical  and  the  3-axis  was  north.  *  The 
input  axes  of  the  gyros  and  accelerometers  were  aligned  along  the  platform 
axes.  The  error  sources  considered  were  initial  position  (3),  initial  platform 
alignment  (3),  accelerometer  bias  (3),  accelerometer  scale  factor  (3)  and 
gyro  bias  drift  (3).  To  evaluate  the  impact  accuracy  of  a  thrust  termination 
at  400  sec,  the  trajectory  tape  was  aborted  ai.  400  sec  and  the  equations  of 
motion  were  used  to  integrate  during  free  flight;  they  were  terminated  when 
the  attitude  west  through  aero  on  a  negative  slope.  Although  the  trajectory 
tape  had  only  one  file,  it  was  necessary  to  enter  a  1  in  TRAJ  in  order  to 
obtain  multiple  processing  of  the  trajectory.  The  data  sheets  used  for  this 
run  are  shown  as  B-l,  Error  Source  Schedule  and  B-2,  Orientation  and 
Control  Data  in  Appendix  B.  Since  two  more  cases  were  to  be  run  by 
ERAN  before  being  processed  by  0UTP,  the  ENDJ0B  O  card  was  scratched 
out  in  B-2. 

The  processing  of  this  data  was  controlled  by  0UTP  and  the  data  sheet 
used  is  shown  as  B-9.  Since  it  was  desired  to  process  the  ERAN  tape  in 
sequential  order,  it  was  not  necessary  to  enter  anything  in  CASE.  If  process¬ 
ing  in  a  different  order,  or  reprocessing  any  narticular  case,  had  been 
wanted,  it  could  have  been  done  by  using  the  CASE  control.  It  was  required 
to  obtain  output  in  ECI  coordinates  at  the  powered  flight  termination  and  at 

—————— - 

Since  the  trajectory  was  run  on  a  spherical  earth  model,  the  geodetic  and 

geocentric  latitudes  are  equal. 


impact,  and  to  evaluate  the  impact  errors  bv  using  mission  evaluation 
Option  1.  As  it  was  desired  to  save  the  SIGMA  data  for  the  next  case,  the 
S0PT  option  was  called  for.  The  l<r  errors  for  this  case  are  shown  in 
Table  3. 

Table  3.  One-Sigma  Errors  for  Test  Case  1 


Vector  Number(e) 

Error  Source 

Sigma  Value 

1,  2,  3 

initial  position 

500  ft  (three  a xea) 

4,  5,  6 

initial  platform  orientation 

30  s'ec  (three  axes) 

7.  8,  9 

accelerometer  bias 

10-4g  (3  components) 

10,  11,  12 

accelerometer  scc’e  factor 

10~4g/g  (3  components) 

13,  14,  25 

gyro  bias  drift 

0.  1  dcg/hr  (3  components) 

The  bias  and  scale-factor  error  sources  of  each  accelerometer  were  corre¬ 
lated  with  a  correlation  coefficient  of  0.  5,  i.  e.  ,  the  number  one  accelerom¬ 
eter  bias  error  (7)  was  correlated  with  its  scale-factor  error  (10),  etc.  Since 
additional  cases  were  to  be  processed  by  0UTP,  the  ENDJ0B  O  card  was 
scratched,  completing  the  input  for  thia  case. 


4.2 


TEST  CASE  2 


In  this  case,  an  evaluation  is  made  of  the  altitude,  cross-range,  and  time 
errors  at  a  fixed  range  after  completion  of  one  orbit.  The  configuration  of 
the  inertial  navigatox  was  the  same  as  in  Ttsst  Case  1,  but  it  was  arr./ed  at 
in  a  different  way.  The  platform  1-axis  was  aligned  with  the  vertical  as 
before,  but  the  3-axis  was  aligned  east.  To  retain  the  same  gyrb  orientation 
with  respect  to  the  trajectory,  it  was  necessary  to  rotate  the  gyro  cluster 
90*  about  its  1-axis.  The  accelerometer  alignment  was  controlled  by  using 
the  nonorthogon&l  accelerometer  option  (see  B-4).  The  error  sources 
considered  were  the  same  as  in  Case  i,  with  the  addition  of  terminal  errors 
(applied  at  the  abort  time).  The  trajectory  tape  was  aborted  at  500  sec  to 
insert  the  terminal  condition  errors  (thrust  tailofi,  guidance  equations,  etc.  ) 
and  the  equations  of  motion  used  for  one  orbit  (approximately  5500  sec).  It 
was  desired  to  have  a  time  history  output;  therefore,  0UT  was  made  non  zero, 
PPF  was  set  at  400  and  PFF  at  2000.  The  data  sheets  used  to  make  this  run 
are  shown  as  B-3  through  B-5  in  Appendix  B.  It  was  also  necessary  to 
scratch  the  ENDJ0B  card- 

The  data  sheet  for  output  processing  of  this  case  is  shown  as  B-10. 

Using  Option  3,  t  me  history,  results  in  additional  output  at  400  sec  during 
powered  flight  and  every  2000  sec  during  orbit.  It  was  desired  tc  have  only 
the  LH  coordinate  system  output;  therefore,  six  zeros  (000000)  were  entered 
in  SYSTM  to  negate  the  logic  from  thj,  Case  l  option.  The  mission  evaluation 
Option  Z  was  used  and  the  S0PT  option  cancelled.  Since  the  SIGMA's  for  the 
first  15  error  sources  were  held  over  from  Case  1,  only  the  terminal  condi¬ 
tion  error-source  sigma  values  were  required.  Since  no  changes  in  the 
correlation  coefficients  from  Case  1  ’"“re  desired,  entries  in  RH0  were 
not  necessary.  The  ENDJ0B  card  was  scratched,  completing  the  input  data 
required  for  this  case. 


4.3 


TEST  CASE  3 


r 


i 


* 


I 


T'iir  was  an  evaluation  of  the  errors  at  injection  into  the  final  orbit.  The 
gyro  drift  was  saumed  to  have  an  exponential  autocorrelation  function  with 
a  time  constan  A  2  hr.  The  accelerometer  bias  during  the  trajectory's 
final  powere>  equence  had  a  standard  deviation  three  times  larger  than, 
and  uncor related  with,  that  of  the  initial  trajectory  sequences.  The  platform 
l-axis  was  aligned  with  the  geocentric  vertical  and  the  3-axis  was  north.  The 
initial  condition  position  errors  were  referenced  to  platform  axes.  The 
platform  was  torqued  about  its  3-axis  (approximately  missile  pitch  avis)  to 
maintain  a  small  (<20° )  angle  between  the  missile  and  the  platform  axis 
throughout  the  trajectory.  The  accelerometer  and  gyro  alignment  with 
respect  to  platform  axes  was  the  same  as  in  Case  1.  The  error  sources 
were  the  eame,  with  the  addition  of  a  gyro-torquing  scale-factor  error. 


To  achieve  the  evaluation  described  above,  the  trajectory  was  divided  into  3 
phases,  with  the  first  terminating  in  2  hr,  the  second  in  4  hr,  and  the  third 
phase  at  the  end  of  the  trajectory  (5.  63  hr).  This  allowed  an  approximation 
of  the  effect  of  the  gyro -error  autocorrelation  function.  (More  phases  would 
more  nearly  approach  the  true  effect.  ) 


The  changes  in  the  Error  Source  Schedule  were  in  EAOO,  EGOO,  EG06,  and 
the  terminal  condition  errors  (see  B-6  in  Appendix  B).  The  body  turning 
rites  are  given  in  B-7,  and  the  changes  in  the  control  data  (B-8)  were 
the  following: 


PSIP 

PHIP 

PSII 

I  GYRO 
I  ACCEL 


0,  aligns  platform  3-axis  north  (4;  =  0) 

-0,  changes  platform  alignment  option  back  to  1  and  aligns  u 
with  respect  to  geocentric  vertical 

-0,  changes  IC  option  back  to  1  and  causes  initial  position 
errors  to  be  along  platform  axes 

0,  references  gyros  to  platform  axes 

0,  references  accelerometers  to  platform  axes 
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I  10  =  0,  changes  accelerometer  option  back  to  1,  i.  e.  ,  orthogonal 

orientation 

TSUBA  =  chosen  to  be  a  time  greater  than  (or  could  have  been  equal  to) 
the  end  of  the  trajectory  tape 

ENOC  =  0,  eliminates  use  of  the  equations  of  motion 

0UT  =  0,  eliminates  use  of  the  intermediate  tape  writing  intervals 

BhST  =  1,  indicatee  a  non-inertial  platform  case 

BRTAB  =  1,  indicates  rates  are  supplied  by  an  input  table 

TG0P  =  7200,  indicates  time  to  end  the  first  phase 

1  =  14400,  indicates  time  to  end  the  second  phase 

The  final  phase  is  ended  by  the  termination  control. 

The  data  sheet  for  the  output  processing  of  this  case  is  given  in  B-ll. 

The  print  option  was  2  (phase  and  trajectory  discontinuities)  and  EVALU  set 
at  zero  to  eliminate  its  output.  Although  the  SIGMA's  were  the  same  for  the 
first  15,  they  had  to  be  re-entered  because  the  S0PT  option  was  zeroed  out 
in  the  previous  case.  The  vector  errors  for  this  case  are  si  own  in  Table  4. 

Table  4.  Vector  Errors  for  Test  Case  3 


Vector  Numbers 

Error  Source 

Sigma  Value 

1,  2,  3 

initial  position 

500  ft  (3  axes) 

4,  5,  6 

initial  platform  orientation 

30  s'ec  (3  axes) 

7,  8,  9  (phase  i) 

accelerometer  bias 

10-4g  (3  components) 

10,  11,  12 

accelerometer  scale  factor 

10'^g/g  (3  components) 

13.  14,  15  (phase  1) 

gyro  drift 

0.  1  deg/hr  (3  components) 

16,  1? 

gyro  torquer  scale  factor 

0  (Nos.  1  and  2  gyros) 

18 

gyro  torquer  scale  factor 

10 “4  (No.  3  gyro) 

19,  20,  21  (phase  2) 

gyro  drift 

0.  1  deg/hr 

22,  23,  24  (phase  3) 

accelerometer  bias 

3X10“4g  (3  components) 

25,  26,  27  (phase  3) 

gyro  drift 

0.  1  deg/hr 

28,  29,  30 

terminal  velocity  errors 

0.  1  ft/sec 

The  correlation  of  accelerometer  bias  and  scale  factor  during  Phase  l  was 
assumed  to  He  the  same  as  in  Cases  1  and  2,  but  it  had  to  be  reentered 
because  additional  correlation  coefficients  were  being  entered.  The  tim%. 
correlation'  of  gyro  errors  were  calculated  as 

-  tj) 

Pij s  exp  -  -jznr 


where  p..  is  the  appropriate  correlation  coefficient,  and 

"  J 

(t.  -  t.)*  =  14400  -  7200  =  V*C0 

l  j' 

20288  -  7200  -  13088 
20288  -  14400  =  5888 


This  completes  the  description  of  the  input  data  required  for  this  case. 

Printouts  listing  the  cards  used  for  ERAN  data  ar«i  0UTP  data  are  given  for 
each  run.  Along  with  the  output  listings,  they  are  included  in  Appendix  C. 


Finer  or  coarser  time  intervals  could  have  been  chosen  to  approximate  the 
autocorrelation  function  exp  -  with  more  or  less  phases  required  for  the 
approximation. 
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Table  A- 5.  ERAN  Equatic;  .  of  Motion  Initialization 
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.  0UTP  Case  Control  and  Data  Input  Form 
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Table  A-7.  0UTP  EVALU  3  Input  Data  Form 
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If  the  roll  and/or  yaw  rates  are  all  sero. 
It  1>  not  necessary  to  enter  N  sero*  for 
the  roll  and/or  /aw  rates.  Slaply  enter  a 
Z  in  the  prefix  field  and  the  mseber  "W"  in 
the  value  field  to  genur  ,te  N  sero#. 


B-7.  Turning  Rates  in  Platform  Coordinates  for  Test  Case  3 
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R-9.  Case  Control  and  Data  Input  Form  ior  i  es'  Case  1 
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APPENDIX  C 

OUTPUT  LISTINGS  FOR  SAMPLE  CASES 
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C-l.  ERAN  Input  Data  ....  C-3 


C-2.  Test  Case  1 .  C-5 
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C-2.  Test  Case  2  (Continued) 
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Test  Case  1  (Continued) 
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C-3.  Test  Case  2  (Concluded) 
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C-4.  Test  Ca3e  3  (Continued) 
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C-4.  Test  Case  3  (Continued) 
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C-4.  Test  Case  3  (Continued) 


phase  logic  syn  eo  mission  TRAJECIORY  TAPE  5*2 


0»M040  •  «■*  •  •* 


OOONSN  l<\  *  (*»  P* 


o 


e  —  r»  <o  f*  Y' 

00-#»»00000000<000 
I  ir^i  til  ii 

X 

o 

« 


oo««mo  00*0000  —  -POO 
*#•##•••••••*••• 

0001*000  —  ©  O  ©  O  —  NOO 
4  11  till  ' 


(O^JOOO^NO^iAOIMOO'-'/' 

«#•#•••*•••••••• 

ooooomoo-^-^ooooro^ 

>  I  I 

> 

o 


x 

UJ 

*- 

Si 

> 

Si 


000000^0’0^0^000‘T> 
•  «•••••••••••••• 

OOOOON^OONOOOOH-4 

K  I  I  I  I  I  I  II  I  I  I  I  I 

> 

3 


h«n01AO''0^^ 
»  w  o  m  O'  i  30 
V  *  OD  -'!’  ^ 

I  I  I 


TN4  m 

f*  0-4 

I 


t*\  <X>  f\»  UTN 

O'  o 
^  in 
1 


K  J  MOIT  I  IT  o  <© 

,0  r»  rg  rg  >f  m  o 


r->f®>f'.n  P-*\J 

r-  -*  *•*  |  — »  rg  in  i  m  ® 

0+  \T<  I  t  ^  <\i  (\f  >f 


o  ®  *1  <NJ  rg 

>$*  I  •  O'  ® 

I  fNJ  nf 

'  i 


oooooooooooooooo 
I  I  I  I  I  I  I  I  •  I  I  I  I  *  I  I 
^rgro-^r^r^^pg^-rfcg^N^cgrifn 
^•^-4r^mpoooo*^^^ooo^ 
M  M  •»  M  —  OOOOOOC  u  O  w 
O  O  O  O  O  ©  * 


o 


C-  52 


Test  Case  3  (Continued) 
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C-4.  Test  Case  3  (Continued) 
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Test  Case  3  (Continued) 
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APPENDIX  D 


RESETS 

0.1  INTRODUCTION 

"Reset*”  is  the  term  applied  to  correction  schemes  for  updating  inertial 
system  navigation  data  as  a  result  of  measurements  made  by  sensors  other 
than  the  inertial  elements  (gyros  end  accelerometers).  These  measurements 
can  be  made  by  ground  equipment  from  which  processed  navigation  data  is 
transmitted  to  the  navigation  computer;  or  airborne  sensor  data  can  be  used 
by  the  computer  to  update  the  navigation  data.  Many  sensor  configurations 
and  data  processing  schemes  could  be  useo,  but  only  a  few  are  discussed 
here.  It  is  assumed  for  this  analysis  that  only  a  limited  number  of  corrections 
would  be  made,  although  many  measurements  might  be  used  and  pr*-filtered 
to  obtain  the  measured  value  used  for  reset. 


PRECEDING 
PAGE  BLANK 
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MEASUREMENTS  CONSIDERED 


D.  2 

The  measurements  considered  are  component*  of  position,  velocity,  and/or 
platform  orientation  as  follows: 

a.  Position  Measurements 


(!) 

altitude 

(2) 

slant  range  to  a  ground  station 

(3) 

position  vector  in  the  ECI  or  orbit  plane  (see  Section 
2,  4,  1)  coordinate  system  as  derived  by  a  ground 
station  and  transmitted  to  the  airborne  computer 

b. 

Veloc 

ity  Measurements 

(1) 

altitude  rate 

(2) 

slant  range  rate  to  a  ground  station 

(3) 

velocity  vector  (ECI  or  orbit  plane  coordinates)  as 
derived  by  a  ground  station 

c. 

Angular  Measurements 

(1)  angular  measurements  with  re  >pect  to  inertial  space, 
using  a  steU«.r  sensor 

(2)  angular  measurements  with  respect  to  the  local 
vertical  (assumed  geocentric),  using  a  horizon  sensor 

(3)  piatlorm  error  vector  as  derived  from  multiple 
stellar  sensor  or  horizon  sensor  measurements. 

In  developing  the  appropriate  reset  equations,  it  was  desired  to  retain  the 
common  notation  of  Y  for  measurement  vectors  and  X  for  state  vectors. 

For  this  reason,  there  is  some  ambiguity  between  this  otation  and  that,  of 
coordinates  of  the  position  and  velocity  vectors  (e.  g.  ,  see  Section  D.  6).  It 
is  hoped  that  this  am'oiguity  will  not  cause  confusion  in  interpretation. 

Three  reset  methods  are  considered  that  use  the  same  form  for  correcting 
the  state  vector.  The  concept  for  deriving  the  form  of  t  ,e  reset  equations  is 
obtained  from  the  following  considerations.  A  rr.easu'ement  is  functionally 
related  to  the  navigation  d  ta  and  is  corrupted  by  some  noise.  This  relation¬ 
ship  is  shown  by 

Y.  =  .F .  ( X )  +  N. 

l  '  l 
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wh  ere 


v  •  th 

Yi  is  i  measurement 

Fj(X)  is  the  functional  relationship  (see  Section  D.  6  for 

explicit  forms) 

X  is  the  true  state  vector  (9  elements) 

th 

Nh  is  the  error  (noise  plus  bias)  of  the  i  measurement 
(a  random  variable) 

Likewise  a  vector  of  measurements  can  be  constructed  as 

Y  =  F(X)  +  N 

The  measurement  state  vector  can  be  estimated  from  the  navigation  system 
data  as 


Y  =  F(X) 

where  X  is  the  navigation  system  estimate  of  the  state  vector  X  and  is  in 
error  by  AX,  that  is 


X  =  X  +  AX 


where  AX  is  a  random  variable,  which  is  the  sum  of  the  effects  up  to  time 
(t)  of  all  the  error  sources  (P  in  number),  that  is 


Then  the  difference  between  Y  and  Y  can  be  calculated  and  used  to  estimate 
AX  and/or  That  is 


AY  =  Y  -  Y 

=  F(X  +  AX)  -  F(X)  -  N 

F(X  +  AX)  -  F(X)  can  be  expanded  into  a  Taylor  Series  so  that 

3F 

AY  =  AX  +  •  -  •  -  N 


where  9F/9X  is  the  matrix  formed  by  taking  the  partial  derivatives  of  F(X) 
with  respect  to  X  for  each  measurement.  All  higher-order  terms  are 
assumed  to  be  zero. 

Again,  for  purposes  of  error  analysis,  it  i6  desired  to  obtain  the  sensitivities 
of  AY  to  each  error  source.  Therefore 

y  =  Mx  -  u 


where  the  notation  has  been  adopted  that 

y  =  sensitivity  (vector  for  multiple  measurements)  of  the 
measurement  state  with  respect  to  the  ith  error  source 
8AY/£€. 

M  =  9F/8X  is  an  m  x  n  matrix  in  which  m  equals  the  number 
of  functionally  independent  measurements  at  time  t  and 
n  =  9  '  — 

x  =  navigation  sensitivity  state  vector  3X/9«. 

u  =  unit  vector  (one  for  each  measurement),  as  described 
in  Section  2.  3.  3.  (An  explicit  symbol  EKlm  is  not  given 
here. ) 


I 
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The  Reset  Equation  is  de/eioped  by  letting 


XR  £  X  -  KAY  (or)  AXR  =  AX  -  KAY 

where  Xp  is  the  estimated  state  vector  after  the  measurement  correction(s) 
and  K  is  an  n  x  m  matrix  developed  in  the  sequel. 

Taking  ;he  partial  derivatives  of  this  expression  with  respect  to  each  error 
source,  it  is  seen  that  the  state  vector  sensitivity  for  each  error  source  is 
changed  when  a  measurement  is  made,  that  is 

6AXR  =  6aX  -  K5AY 

or  xR  =  x  -  Ky 

=  x  -  KMx  +  Ku 
=  (I  -  KM)  x  +  Ku 

NOTE:  The  notation  and  form  of  this  equation  are  the  same  as  for  Kalman 

filtering  (Reference  6),  except  that  x  and  u  are  f  cnsitivity  (w.  r.  s„  t. 
random  variable)  vectors,  rather  than  vectors  of  random  variables. 

It  remains  to  establish  the  matrix  K  in  the  above  equation.  Three  me  the  Is 
are  considered  here,  which  are  termed:  replacement  rose'.,  deterministic 
reset,  and  linear  statistical  reset. 

The  development  of  each  of  these  metnods  will  now  be  descr;bed. 
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REPLACEMENT  RESET 


D.  3 


In  using  this  technique,  it  ia  assumed  that  the  measurement  error  is  much 
smaller  than  the  state  vector  error;  thus,  the  measurement  difference(s) 
(AY)  is  assumed  to  be  dependent  on  the  navigation  state  vector  error  only. 
The  constraint  of  the  correction  is  such  as  to  make  AYn  =  0,  and  therefore 


AYr  =  Yr  -  Y  =  0  =  MAXr  -  N 


Taking  the  partial  derivatives  of  this  expression  results  in 

Vr  =  Mxr  -  " 

=  (I  -  MK)(Mx  -  u)£  0 

In  order  for  this  cond’ti'-'"  to  be  satisfied  for  all  x,  the  first  term  must  equal 

zero.  Therefore  _ 

Li 

MK  =  I 


X 

post  multipling  by  [MM  J 

MK[  MM 1  j  =  mmt 

from  which  K  is  determined  as 

Tr  T  -  1 
K  -  M  [  MM  i 

K  is  sometimes  referred  to  as  the  pseudo-inverse  of  M. 

Thus,  the  reset  sensitivity  state  vectors  for  all  error  sources  active  up  to 
time  t,  based  on  this  constraint,  are 

xr  =  (I  -  KM] x 
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Additionally,  new  sensitivity  state  vectors  are  generated  tc  account  for  each 
measurement  error  by 


x  *  Ku 
m 

which  are  essentially  initial  condition  errors  at  time  t.  and  x  represent 

*  m 

the  total  set  of  sensitivity  vectors  x,  which  are  subsequently  operated  upon  as 
independent  sensitivity  vectors.  Thus,  any  operations  including  additional 
resets  at  time  t  would  include  xm  in  the  set  of  state  vectors. 

X  *1 

A  few  remarks  about  [MM  ]  are  in  order  here.  For  most  measurement 
types  considered  (e.g.  ,  an  individual  po  i»iion  measurement,  a  stellar  '“ensor 
sighting  or  a  platform  reset,  or  a  position  and/or  velocity  correction  from 
the  ground),  [.MM ^  ’  *  =  I. 

For  individual  velocity  measurements  and  nonorthogonal  measurements  (e.  g.  , 
simultaneous  altitude  ard  ulant  range  measurements),  [ MM  ]  I. 

Functionally  dependent  measurements  (such  as  altitude  and  slant  range  when 
the  slant  range  vector  is  along  the  radius  vector,  and  overdetermined 
measurements  when  two  stellar  sightings  give  four  measurements  of  platform 
angle  errors),  result  in  singular  inverses.  The  pseudo-inverse  could  be 
invoked  for  the  altitude  and  slant  range  example,  which  would  in  effect  result 
in  the  average  of  the  two  measurements.  The  case  of  the  twc  stellar  sightings 
should  be  reformulated  by  partitioning  the  state  vector  into  the  orientation 
elements  only,  from  which  the  least  squares  solution  {[M  M]  M  )  can  be 
used. 

The  question  of  sequential  (at  time  t)  reset  vs  simultaneous  reset  is  of 
interest.  The  M  matrix  is  formed  from  row  vectors  of  partials.  If  the 
dot  product  oi  these  vectors  is  zero,  and  the  cross  product  is  one,  then 
sequential  or  simultaneous  'reset  is  equivalent.  However,  if  both  these 
conditions  are  not  satisfied,  the  order  of  sequential  reset  is  important. 
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A h  eU*mt«  and  generally  equivalent  approach  is  one  in  which  a  transformatiu 
matrix  ia  developed,  which  transforms  the  state  vector  in  EC  I  coordinates 
into  a  Stats  vector  in  ths  measurement  coordinate  system.  The  elements 
that  are  measured  are  set  to  aero,  and  the  resultant  state  vector  is  trans¬ 
formed  back  by  the  inverse  of  the  transformation  matrix.  This  method 
suffers  in  that  simultaneous  multiple  measurements  usually  cannot  be  included 
in  the  reset  technique;  therefore,  multiple  measurements  would  nave  to  be 
handled  sequentially  and  sometimes  only  in  a  specific  order. 


:  f 
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D.  4  DETERMINISTIC  RESET 

In  this  method  it  is  assumed  that  the  measurement  differences  are  solely 
dependent  on  an  equal  number  of  error  source*.  A  typical  example  would  be 
the  use  of  a  stellar  tracker  to  derive  the  launch  position  errors  of  a  mobile 
missile  system,  and  thus  compensate  for  the  error.  The  technique  developed 
here  is  not  restricted  to  measurement  types  or  the  error  sources  to  be  con¬ 
sidered.  The  concept  is  developed  as  follows.  Let 


where 

AXjj  =  the  error  in  X  due  to  the  error  source  vector 

D  =  an  r.  x  m  matrix  formed  from  the  sensitivity  vectors 
°f,D 

«  =  an  m  vector  of  error  sources  to  be  determined  from 

the  measurement  s) 

In  using  this  technique,  it  is  assumed  that  the  measurement  error  is  much 
smaller  than  the  effect  of  the  state  vector  error,  and  that  the  state  vector 
error  due  to  all  other  error  sources  is  much  smaller  than  the  state  vector 
error  due  to  i^.  Therefore,  the  measurement  differences  (AY)  are  assumed 
to  be  dependent  on  the  error  vector  only,  so  that 

4Yd  =  MAXd  .  MD.d 
<D=  [MD]-1AYd 
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However,  tD  can  only  be  estimated  from  AY  (the  quantity  derived  from  the 
meaeuremtnte),  so  that  the  estimated  value  of  «D  is 

€d  *  [MD]-1AY 

The  state  vector  estimate  is  corrected  by  using  the  estimate  of  «D, 
as  follows 


xR  =  x  -  d«d  =  x  -  d[md]_1ay 

Thus,  it  is  seen  that  for  this  correction  scheme 

k  =  d[md]_1 


and 

G 

xr  s  (I  -  KM)x 

The  sensitivity  vector  corrections  follow  the  same  procedure  as  for  the 
replacement  reset;  each  sensitivity  vector  is  corrected  at  time  t,  and  new 
sensitivity  vectors  are  added  to  account  for  the  measurement  errors.  The 
constraints  on  [MD]  *  are  that  the  measurements  be  functionally  independent 
and  that  the  error  sources  to  be  determined  do  not  have  equal  sensitivities. 

It  is  apparent  from  the  state  vector  reset  equation  that  the  explicit  determina¬ 
tion  of  is  not  required  for  corrections  of  the  navigation  data.  Also,  it 
can  be  shown  that  the  reset  sensitivity  state  vector  for  each  element  of  the 
vector  is  set  to  zero  at  the  reset  time,  provided  that  the  sensitivity 
matrix  D,  used  in  the  reset  equations,  is  formed  from  the  sensitivity  vectors 
of  the  error  analysis  discussed  in  detail  in  the  next  paragraph.  Thus,  the 
navigation  errors  resulting  from  vector  errors  are  zero  at  reset  time, 

(  '< 


D-  12 


regardless  of  the  magnitude  of  the  vector  (assuming  the  elements  of 
are  not  large  enough  to  invalidate  the  linearity  assumption).  If  the  vector 
were  composed  of  initial  condition  elements  only,  the  effects  would  be  negated 
and  new  error  sources  formed  by  the  measurement  error  vector.  On  the 
other  hand,  if  the  vector  had  elements  representing  forcing  functions 
(such  as  accelerometers,  drag,  etc.),  it  might  be  desirable  to  use  the 
estimate  of  to  compensate  those  parameters,  thereby  reducing  the 
variance  of  navigation  errors  following  the  reset  time.  For  this  procedure, 
it  is  necessary  to  determine  the  variance  of  the  estimate  error,  which  can 
be  calculated,  when  desired,  from  the  relationship 

«D  =  [MD]"lAY 

=  [MD]_1{MAXD  +  MAXp_D  -  N) 

P-D  M 

3  «p  +  [  MD]  ~  *  M  6  y;«i  -  [MDl^yy 

i=l  j=l 


where 

are  all  error  sources  excluding  the  set 

«.  are  the  measurement  errors 
J 

Thus,  the  estimate  error  is  defined  as 

AtD  =  *D  ‘  *D 

and,  assuming  that  the  measurement  errors  are  independent  of  the  other 
error  sources,  the  variance  covariance  of  the  estimate  error  can  be 
calculated  as 

e<A<dV)  =  +  q|([MD]*ST 
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where 

X  is  the  covariance  matrix  of  navigation  state  vector  errors 
'  excluding  «_  parameters  at  reset  time  (before  the  mea¬ 
surement  ir  included) 

Q  is  the  covariance  matrix  of  measurement  errors  at 
reset  time 

In  the  application  of  this  method,  the  sensitivity  matrix  D  would  be  obtained 
in  one  of  the  three  following  ways: 

a.  Input  based  on  a  nominal  trajectory 

b.  Computed  by  using  the  normalized  integral  approach 
(which  excludes  the  effects  of  gravity  feedback  in 
navigation  equations) 

c.  Computed  in  a  manner  similar  to  the  equations  of  this 
program. 

Appi  ^ach  (a)  should  be  adequate  for  most  application  i  for  the  same  reason 
that  a  nominal  trajectory  suffices  for  error  analysis.  Approach  (b)  more 
closely  approximates  the  true  sensitivities  for  a  given  trajectory  and  is  not 
a  difficult  computation  for  an  airborne  computer.  An  analysis  of  either  of  the 
first  two  methods  requires  the  generation  of  the  D  matrix  from  alternate 
runs,  which  would  then  be  treated  as  input  data.  Method  (c)  is  self-contained, 
but  more  complex  for  airborne  computations. 
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D.  5 


UNEAR  STATISTICAL  RESET 


0 


This  is  equivalent  to  the  Kalman  filtering  technique  (Reference  6).  but  it  is 
developed  differently.  Certain  assumptions  are  made  to  facilitate  the  develop¬ 
ment  of  the  problem  at  hand*  but  as  the  technique  presented  is  a  general  one, 
assumptions  need  not  be  made.  The  one  basic  restriction  made,  on  both  the 
Kalman  filtering  technique  and  the  method  here,  is  that  of  the  best  linear  esti¬ 
mate  of  the  random  variables  in  the  sense  of  minimizing  the  mean  square 
error  of  the  estimate.  The  development  uses  the  conditional  expectation  func¬ 
tion  directly.  For  the  scalar  case  it  is  well  known  that 


t 


a 


e(x|y) 


(The  notation  E  is  adopted  to  distinguish  it  from  the  normal  operator  E  in 
the  conditional  expectation  function,  because  of  the  restriction  of  best  linear 
estimate.  However,  if  the  stochastic  processes  are  Gaussian,  the  operators 
are  equivalent.  This  same  noi&tion  is  used  in  Reference  6. ) 


£(x|y)  = 


E(XY)<rv  , 

- —  *  y  =  E(XY)[  £(Y4)] 

<rx<rY 


-i 


Y 


For  X  and  Y  vectors,  it  can  be  shown  that  the  same  form  holds  as 


E(X  |  Y)  =  E(XYT)[  F(  YYT)]  "  1  Y  =  KY 

where  the  elements  of  K  in  this  expression  (an  n  x  m  matrix)  are  termed 
regression  coefficients  (Reference  7,  Chapter  23,  p  302). 


C 
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For  the  problem  at  hand,  it  is  desired  to  estimate  AX,  the  error  in  the  state 
vector  estimate,  given  the  measurement  difference  AY,  or 

AX  *  E(AX  |AY)  *  E(AXAYT)l  E(AYAYT)  J  ' 1  AY 

and  based  on  this  estimate,  the  state  vector  estimate  would  be  corrected  as 
follows 

XR  *  X  -  AX  =  X  -  KAY 

where 

t:  =  E{AXAYT)lE(AYAYT)]"1 

As  defined  previously 


AY  =  MAX  -  N  ,  AYT  *  AXTMT  -  NT 


and  therefore 

K  =  E(AXAXTMT  -  AXNT)[M  E(AXAXT)MT  -  M  E(AXNT)  -  E(NAXT)MT 

+  E(NNT)]"J 

Assuming  that  the  inertial  navigation  system  parameters  are  independent  of 
measurement  errors,  this  reduces  to 

K  =  J$mT  '  DnRn(t’  t)([M^mT  '  M  DnRn(t'  '  Rl{t>  t)DJmT  +  Rn(t>t)]’1 


D-16 


where 

^  is  the  covariance  matrix  of  navigation  error*  at  reeet  time 

R  (t,  t)  is  a  matrix  of  time  correlation  functions  {including 
cross-cor  relation  terms)  for  the  measurement  errors 

Dn  is  a  matrix  that  propagates  the  effects  of  measurements 
made  at  a  previous  time  (or  at  the  same  time  but  processed 
sequentis 1  ly) 

Rn(t,  t)  represents  the  covariarce  matrix  of  measurement 
errors  at  time  t. 

If  it  is  assumed  that  the  measurements  are  independent  (and  time  and  time- 
cross-correlation  functions  are  zero),  the  gain  expression  can  be  further 
reduced  to 


k  =  ^MT[m^MT  +  q]'1 


The  effects  of  the  previous  assumptions  can  be  determined,  however,  since 
the  calculation  of  ^  in  the  program  includes  the  effects  of  time  and  time- 
cross- cor  related  errors  independent  of  the  assumptions  for  obtaining  K,  The 
simplified  gain  computation  represents  a  more  realistic  one  for  airborne  use. 
Finally,  the  reset  state  vector  sensitivity  is 

xr  =  x  -  Ky 
;=  [l  -  KM]  X 

which  is  in  the  same  form  as  the  previous  methods.  Additionally,  a  new 
vector  is  formed  for  eacn  measurement,  and  the  complete  set  of  sensitivity 
vectors  is  handled  in  the  same  manner  as  Methods  1  and  2  (in  Sections  D-3 
and  D-4). 
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This  reset  technique  could  b«  extended  to  include  the  estimation  of  the  es  /or 
eourcee  a*  follow* 

W«  |AT)  *  Et«AYT)lE(AYAYT)]  “  lAY 

where  *  *  the  error  vector  to  be  estimated,  and 

AY  *  MAX  -  N 

»  MD*  -  N  ■  H«  •  N 

where  D  =  the  sensitivity  matrix  of  all  error  sources. 

Therefore 

«  =  E(«  IaY)  =  e[«  «tht  -  «NT][E(AYAYT)]"lAY 

T 

Making  the  same  assumptions  as  before,  that  E(«N  )  =  0, 

«  =  ^ht[m^mt  +  q]-1Ay 


where 


=  E(*  «7)  *  covariance  matrix  of  all  error  sources 
can  also  be  written  as  D 


By  using  the  above  relationships,  c  can  also  be  written 


*  =  [  DTD]  " 1  DTKAy  =  Kt  Ay 


T  - 1 

provided  that  [D  O]  exists.  If  it  is  singular,  then  the  pseudo* inverse, 
which  always  exists,  can  be  developed.  However,  it  is  not  required,  since 
the  original  formula  does  not  present  this  problem. 
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For  purpose*  of  an  error  analysis,  it  is  necessary  to  derive  the  variance  of 
the  estimate  error  =  c  -  «),  which  is 


=  E{(c  -  «)<«T  -  *T)} 


and  reduces  to 

where  i*  the  covariance  matrix  of  the  reset  error  source  parameters. 

There  are  several  approaches  that  can  be  used  to  compensate  the  navigation 
data,  given  the  estimate  of  the  error  vector;  the  most  straightforward  being 
to  restfi  the  sensor  error  equation. 


It  is  seen  that,  in  theory,  the  navigation  data  and  error  source  parameters 
can  be  corrected,  given  one  or  more  measurements.  Thus,  from  the  stand¬ 
point  of  minimum  navigation  error  in  the  sense  of  the  least-mean- squared- 
error  under  the  constraint  of  a  linear  estimate,  an  optimum  use  of  the  data 
would  incorporate  both  navigation  data  and  error- source  :  ^set. 

Practically,  the  task  described  would  overburden  any  airborne  computer,  and 
in  particular  the  computation  of  the  D  matrix.  Although  the  estimate  of  t 
could  be  partitioned  so  that  only  a  selected  few  would  need  to  be  estimated, 
all  sensitivities  are  required  {or  should  be)  for  the  calculation  of  For 
that  reasc  ■  it  is  assumed  that  calculation  of  the  sensitivities  (D  matrix 
elements)  would  be  based  on  the  nominal  trajectory;  therefore,  it  would  be 
precomputed  and  inserted  in  the  airborne  computer  for  flight  data  processing. 

In  the  special  case  of  free  flight  (orbit  navigation),  most  forcing  functions 
reduce  to  zero;  and  the  problem  can  be  reformulated  do  that  ]j£  can  be  updated 
by  using  the  transition  matrix  in  conjunction  with  those  forcing  functions 
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acting  in  orbit  (drag,  gravity  model  constants,  gyro  drift,  control  system 
impulses,  etc.).  The  capability  to  perform  this  last  analysis,  that  is,  orbit 
reset  or  orbit  navigation,  is  presently  under  development  in  a  separate 
program.  A  report  describing  it  will  be  published  by  J.  Meditch  of 
Aerospace  Corporation. 
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M  MATRIX  GENERATION 


D.  6 

The  M  matrix  is  assumed  to  be  computed  by  the  airborne  computer  and  based 
or  the  functional  relationships  of  the  measured  quantities  (Y)*  and  the  naviga¬ 
tion  data  (X,  X,  Mgp).  Each  measurement  represents  a  row  of  the  M  matrix 
and  is  developed  as  follows  for  the  measurements  considered. 

D.  6.  1  Altitude 


Y  =  h  =  R  -  R  =  JX1  +  Y2  +  Z2  -  R 


6Y  -  f>X  =  65T 


and 


m 


X  Y  Z 

0  5 


0  0  0  0  0 


0 

- 


where  it  has  been  assumed  that  Bk^/dX  =  0 
D.  6.  2  Slant  Range  to  a  Ground  Station 

The  slant  range  vector  is  defined  as 


S=Kg-K 

=  <xG  -  x>Xu  + :  <*G  '  Y)7u  +  {zc  -  z)7u 


xs7u  +  YSYU  +  ZST 


u 


5c  ,  # 

See  comment  on  notation  at  end  of  Sec  Hen  D-£. 
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where 


XG  -  \  coe  \Q  coa  0Q 
Yg  =  Re  coa  XQ  ain  0Q 


ZG  =  Re  ain  \Q 


R  = 


A(1  -  e) 


J  I  +  (e^  -  2e)  cos“\.. 


and  the  terma  in  thia  equation  are  aa  defined  in  Section  2.  3.  3. 


Y  =  S  =  ./Xg  +  y|  +  Zg 


6Y  =  ||  6X  =  -  6X 


m  = 


X„  Yn  Z„ 
S  S  S 

'  "S  §~ 


000  000 


It  is  assumed  that  station  locations  are  perfectly  known,  and  also  that  timing 
errors  are  zero. 

D.  6.  3  Position  Vector 

D.  6.  3.  1  ECI  Coordinates 

m  =  [I(3x3)°(3x3)°(3x3)] 
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D.  6.  3.  2  Local  Coordinatea 


where  p, 
D.  6.4 


from  which 


m 


=  [m 


LEC{3x3)°(3x3)J 


ia  defined  aa  in  Section  2.  4. 1 
Altitude  Rate 


„  _  l  _  dR  _  «  _  i?  _  XX  +  YY  +  ZZ 
Y  =  h  =  St=  R  =TTX  - R - 


_  yT6y  t  r6r 


1  •  np  •  'r  'p  • 

=  -V[(R3Tt  -  R5TT).6X  +  RXT63f] 

RZ 


m 


•XR  RY-YR  RZ-ZR  X  Y  Z 

? — ? — r"5' 


o 
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Slant  Range  Rate 


where 


*s^u  +  + 

(XG  -  +  (Yg  -  Y)Ty  +  (ZG 

<-“.*G  '  *>SU  +  <“.XG  '  Y)?U 


and  the  form  of  m  is  the  same  as  for  altitude  rate  with  the  substitutions  of 
5T  =  5fg.  5T  =  Xg,  R  =  S,  and  R  =  S.  Station  location  and  timing  errors  are 
assumed  to  be  zero. 

D.  b.  b  Velocity  Vector 

D.  6.  6.  1  ECI  Coordinates 


m  =  [Vxa^pxSjVxS)] 

D.  6.  6.  2  Local  Coordinates 

m  =  [°{3x3)MLE°(3x3)] 

D.  6.  7  Stellar  Sensor  Measurement 

It  is  assumed  that  the  sensor  is  mounted  on  the  platform  and  can  be  slewed  in 
azimuth  about  the  platform  1-axie,  and  in  elevation  about  the  tracker  2-axis 
to  sight  on  a  prescribed  star.  Figure  D-l  shows  the  tracker  axes  system  with 
respect  to  the  platform  axes.  The  transformation  between  a  vector  ir  platform 
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coordinate*  and  tracker  coordinate*  i* 


• 

cet 

0 

-SEt 

c 

1 

0 

0 

mtp* 

0 

1 

0 

0 

CAt 

-SAt 

S£t 

■ 

0 

CEt 

0 

SAt 

CAt 

where 

A_,  is  the  azimuth  of  the  star  expressed  in  nominal 
platform  coordinates 

Ej  is  the  elevation  of  the  star  expressed  in  nominal 
platform  coordinates 

The  tracker  is  capable  of  measuring  coordinates  (small  angles)  only  about  its 
2-axis  (a  coordinate  along  the  T1  axis),  and  about  its  1-axis  (a  coordinate 
along  the  T2  axis).  Alternately,  it  can  measure  the  changes  in  AT  and 
which  make  the  above  coordinates  zero.  Mathematically  either  measurement 
type  is  equivalent.  Since  the  platform  rotation  errors  are  assumed  to  be 
small,  they  can  be  treated  as  vectors,  so  that 


- 1 

H 

_ 1 

CEt 

”  SEjSAj 

-SEtCAt' 

L02T. 

II 

o 

_ 1 

CAX 

-SAt 

where 

0^  is  a  measurement  along  the  sensor's  2-axi» 
02T  is  a  measurement  along  the  sensor's  1-axis 


Therefore,  for  a  single  stellar  fix  with  a  stellar  sensor,  two  angles  can  be 
used  to  correct  the  effec*  of  platform  errors.  In  this  case 


0 

0 


oooo  cet  -setsat  -setcat 

0  0  0  0  o  cat  -sat 


For  the  case  of  two  independent  stellar  sightings  {two  different  stars, 
preferably  90*  apart  and  in  the  platform's  2^3^  plane),  the  platform  orientation 
errors  are  overdetermined.  In  this  case,  least  squares  or  other  techniques 
could  be  U3wd  to  process  the  four  measurements  so  that  the  platform  error 
angles  were  determined  explicitly.  This  is  implied  in  the  last  measurement 
category  (Section  D.  6.  9),  Platform  Error  Vector. 


D.  6.  8  Horizon  Senso  Measurements 

A  horizon  sensor  primarily  measures  small-angle  deviations  of  its  mount 
with  respect  to  local  vertical;  the  measurements  can  also  be  processed  to 
indicate  altitude.  Conventionally,  it  is  mounted  on  the  vehicle  frame  and 
used  as  a  reference  for  the  vehicle  control  system,  while  maintaining  small- 
angle  deviations  of  the  vehicle  axes  with  respect  to  local  vertical. 


In  conjunction  with  platform  gimbal-angle  readout  (direction  cosines  in  the 
case  of  strapped-down  inertial  systems)  and  the  navigation  system  position 
data,  the  angles  measured  by  the  sensor  can  be  used  to  correct  the  platform 
angular  errors,  navigation  system  position  errors,  or  both.  For  purposes 
of  a  general  analysis  of  the  use  of  the  horizon  sensor,  the  local  horizontal 
coordinate  system  will  be  used.  It  is  assumed  that  the  sensor  measures 
pitch  (rotations  about  the  Z-axis  of  the  local  horizontal  system)  and  roll 
(rotations  about  the  X-axis).  Transformed  into  this  system,  the  state  vector 
error  is 

=  MLE^ 

where  M,  —  is  as  described  in  Section  2.  4.  1. 

LE 
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Undar  these  assumptions,  the  eenaor  measures 


0 


A®  *  -  +  f  2  "  n0  pitch  measurement 


a Z 

A#  *  -Jr  +  fx  - 


roll  measurement 


where 


AX  *  range  error  in  local  coordinates 

AZ  =  cross-range  error  in  local  coordina.es 

a  platform  error  about  Z  in  local  coordinates 
=  platform  error  about  X  in  local  coordinates 
ng  and  n  =  pitch  and  roll  sensor  errors. 


Thus 


m  » 
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£ 
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M 


LE 


<  ) 
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Platform  Error  Vec.  r 


In  this  case,  it  is  assumed  that  the  platform  error  angles  are  derived  from 
multiple  stellar  sensor  measurements  and  expressed  in  platform  coordinates. 
Therefore 


111  "  [°{3x3)°(3x3)MPEJ 


The  M  matrix  for  any  given  reset  is  then  constructed,  batted  on  any  one  or 
more  of  the  m  matrices  presented  above.  The  combinations,  however, 
must  be  limited  #o  that  elements  of  the  state  vector  are  not  overdetermined 
when  the  Method  1-type  reset  is  used. 


reset  equation  summary 


mvwmm. 
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This  section  summarizes  the  pertinent  equations  presented  for  reset. 
Those  marked  with  an  asterisk  could  be  conveniently  mechanized  in  the 
error  analysis  program. 


X  =  X  +  AX 


navigation  system  estimate  of 
the  state  vector 


AY  =  Y  -  Y  =  MAX  -  N  navigation  system  processed 

measurement  difference 


XR  =  X  -  KAY 
AXr  =  AX  -  KAY 


navigation  system  reset  estimate 
of  the  state  vector  given  the 
measurement(s)  Y 

state  vector  error  following  reset 


*  xr  -  [I  -  KM] x 


*  x  =  Ku 
m 


*  xx 


reset  error-source  sensitivity 
vector(s)  for  all  error  sources 
active  prior  to  current  reset  time 

added  sensitivity  vector(a)  to 
account  for  the  reset  measurement 
error 

measurement  sensitivity  matrix 
(see  Section  D.  6) 


K  is  determined  by  one  of  the  following  3  methods: 

a.  Replacement 

*  k-mt[mmt]'1 

b.  Deterministic 

*  K  =  Dl  MD]  "  1 

where  D  is  formed  from  the  prescribed  state  vector  sensitivities  or  is  input. 

*  =  [MD]  ’  lAY  navigation  system  estimate  of  the 

u  prescribed  error  sources,  given 

the  neasurementfaj  Y. 


At 


D 


*  D  ' 


D 


error  of  the  estimated  error  sources 


( 


i 
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where 


5ptD  =  mt  +  q|([mo]  " *)T 

covariance  matrix  of  the  estimated 
error  sources 


Kris  the  covariance  matrix  of  the  navigation  state  vector  errors 
•'due  to  all  error  sources  excluding  cq  parameters  at  reset  time 
(before  the  measurement  is  included) 

Q  is  the  covariance  matrix  of  measurement  error(s)  at  reset  time 

c.  Linear  Statistical 


K  = 


T 

M  1  +  Q 


-  1 


wh  ?re 


/[.is  the  covariance  matrix  of  the  navigation  state  vector  error 
'due  to  all  error  sources  at  reset  time 

Q  is  the  covariance  matrix  of  measurement  error(s)  at  reset  time 

and  jJ]  and  Q  are  input. 


f  [MD]  1  M 

+  q 

ft 

“7" 

-I 


AY 


navigation  system  estimate  of  the 
error  sources,  given  the  measure¬ 
ment  Y.  (  could  be  partitioned  so 
that  only  «  selected  few  need  be 
explicitly  derived,  e.  g.  ,  e  as  in 
Method  2 


where 


D  is  the  sensitivity  matrix  of  all  (or  partitioned)  error  sources 


is  the  covariance  matrix  of  all  (or  partitioned)  error  sources 
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’fjtBUSWfcw  i  wwn 


where 


=  the  covariance  matrix  of  error  sources  after  reset 

\ 
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DRAG  ERRORS 


APPENDIX  E 


DRAG  ERRORS 


Iii  this  appendix  is  discussed  the  proposed  method  for  estimating  the  effects 
of  atmospheric  drag  errors,  when  the  accelerometers  are  disconnected  dur¬ 
ing  orbital  flight  phases. 

If  there  were  no  forces  experienced  by  the  vehicle  during  coast  periods 
(parking  orbits,  transfer  orbits,  etc.),  the  most  advantageous  way  to  operate 
the  navigation  system  would  be  to  disconnect  the  accelerometers  during 
these  periods  so  that  the  accelerometer  bias  error  would  not  be  integrated. 
For  low-altitude  orbits  (in  the  region  of  100  n  mi),  aerodynamic  drag  force 
is  not  negligible,  but  is  generally  of  the  same  order  of  magnitude  as 
accelerometer  bias.  Consequently,  it  must  be  decided  either  to  measure 
drag  via  the  accelerometer,  or  to  predict  it  by  an  empirical  formula.  The 
decision  would  be  contingent  on  which  method  would  result  in  the  least 
error,  i.  e. ,  on  the  uncertainty  of  accelerometer  bias  vs  the  uncertainty  of 
drag  calculations.  To  fully  answer  this  question,  a  detailed  knowledge  of  the 
configuration  and  flight  time  (function  of  time  of  day,  month,  and  year)  is 
necessary.  For  purposes  of  the  error  analysis,  these  characteristics  are 
generalized  so  as  to  assess  the  relative  importance  of  drag  and  thereby 
determine  if  more  detail  is  required.  Therefore,  the  configuration's 
ballistic  coefficient  and  the  parameters  of  the  atmospheric  density  model  are 
treated  as  random  variables,  with  assumed  means  and  standard  deviations. 

When  utilizing  a  drag  model  for  calculating  the  sensed  acceleration,  the 
equations  of  motion  become 


+  A 


D 


E-3 


where 


Ap  *  vector  of  drag  accelerations 

(assuming  that  drag  acts  along  the 
negative  inertial  velocity  vector) 

where 

Cq  =  drag  coefficient 

S  =  reference  area  (ft2) 
m  =  system  mass  (slugs) 

2 

q  =  dynamic  pressure  (lb/ft  ) 

=  1/2  pV2 

where 

p  =  atmospheric  density  (slug/ft  "*) 

V  =  magnitude  of  inertial  velocity  vector  (assuming  magnitude 
of  inertial  velocity  equals  magnitude  of  relative  velocity) 

*D  '■  -  <!f  *  -  AD  * 


where 

B  =  =  ballistic  coefficient 

V 

gQ  =  reference  gravity  constant  (=  32.  174) 
W  =  vehicle  weight  (lb) 
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The  linearized  differential  equations  are 


The  third  term  is  simply 


B  - 


W 

CDS 


6C, 


6B  =  -  B 


(6  W  =  6S  =  0) 


'D 


Combining  these  results  into  the  linearized  differential  equations  results  in 


—  —  ,  a  yLr  6p{h  )  6C-.  • 

6X  =  M  JX  -  A n{^-~  +  -  I?  +  — iL  +  ~2.)X  -  An6X 

G  D  v2  Mh  R  P(  ho)  CD'  D 


^Mg  +  M^X  +  M^X+Fa 


where 


M  =.‘4*  vfiSSli 

•  3h  ADVl  RV  ' 


ix  p  9h 


m«='adI37-  tl] 


Fa  =  -  (EQOO)AqX 


In  this  particular  case,  EQOO  is  a  scalar  and  equals  1  when  deriving  the 
sensitivity  of  aerodynamic  forces  (i.  e,  ,  due  to  the  rss  value  of  6p(h  )/ 

When  considering  re1  ttively  short  time  durations  {less  than  2  orbits)  in  a 
drag  environment  (a  100- mile  orbit),  these  expressions  can  be  further 
simplified  as 


E-7 


I 


*•  the  coefficient  of  Mq  is  =  1,4X10"®  and 


“fa-?*  V1-10' 


for  B  in  the  range  10  <  B  <  100 


M.  >  >  M. 

G  xx 


M.  =  C 
xx 


b.  the  coefficient  of  M<«  =  A~.  =  4  X  10" ^/B  <  4  x  10"®  and 

xx  D 


6x  >  >  6x 


eo  that  M_  6x  >  >  >  M. .  6x 
G  xx 


M*  •  =  0 

xx 


The  validity  of  these  assumptions  can  always  be  ascertained  by  utilizing 
Reference  9.  which  includes  these  effects  as  well  as  the  first  earth  oblate¬ 
ness  term  ( J^).  The  intent  here  is  only  to  assess  the  relative  importance 
of  atmospheric  effects,  so  that  when  the  results  indicate  that  more  accuracy 
is  required,  Reference  9  will  be  used  for  detailed  studies.  Thus,  the 
equations  for  drag  uncertainty  reduce  to  the  same  differential  equations  as 
those  for  the  other  error  sources,  with  the  forcing  function  of 


Fa  =  -  (EQOO)AdX 


E-8 


where 


B  = 


PV) 

(an  input  constant) 


log  p  =  ^(h)  +  Kg  (an  input  table  plus  a  scaling  constant) 
p  =  e  2.  30258  log  p 


It  remains  to  establish  methods  for  specifying  the  mean  values  of  B  and  p, 
and  some  estimate  of  their  deviations.  With  the  assumptions  that  the  region 
of  interest  is  in  orbital  velocities  at  greater  than  80  miles,  and  the  configura¬ 
tions  of  interest  are  upper  stages  with  payloads  attached,  the  drag  coefficient 
for  zero  angle  of  attack  and  molecular  flow  can  be  approximated  by 


Gjj  =  Cjjj,  +  0.  256 


L 

D 


where  is  the  drag  coefficient  for  the  shape  of  the  payload  section. 

Typical  values  arra 


flat  plate  Cjjp  = 
20*  cone  Cj-jp  = 
15*  cone  Gjjj,  = 
L  = 

D  = 


2.  11 
2.  04 
2.  03 

the  length  of  the  cylindrical  section 
(excluding  cone  sections) 

the  diameter  of  the  cylindrical  section 


Using  this  formulation  of  drag  coefficient  produces  approximately  a  20- 
parcont  uncertainty  in  its  magnitude,  provided  the  assumptions  of  molecular 
flow  and  sors  angle  of  attack  are  maintained.  f  ■ 

The  upper  atmospheric  density  is  affected  by  many  parameters,  by  far  the 
most  by  solar  activity.  In  Reference  10  are  discussed  the  various  factors 
that  influence  the  density  and  they  are  summarized  as  follows: 

a.  Diurnal  f  Day- Night  Effect^.  The  density  varies  with  the 
time  of  aay  f  having  "it  s  minimum  at  night  and  maximum  at 
approximately  2:00  P.  M,  local  standard  time.  The  effect 
is  a  function  of  the  angle  between  the  earth  sun  line  and  the 
radius  vector  and,  therefore,  depends  on  lattitude  as  well 
as  longitude.  The  effect  is  small  (15  percent)  at  low 
altitude  (100  n  mi)  and  increases  with  altitude  to  more  than 
100  percent  at  200  n  mi  (utilizing  Eq.  (10)  in  the  Jacchia 
formula,  given  in  Reference  1 1). 

b.  Solar  Activity  ^1 1-  Year  Cycle).  There  is  still  a  rather 
forge  discrepancy  in  models  for  this  effect,  as  evidenced 
by  the  curves  presented  in  Reference  10:  Figure  1  for  the 
Jacchia  I960  model  and  Figure  2  for  the  Paetzold  1962 
model.  However,  there  is  reasonably  gooa  agreement  in 
the  low- altitude  region.  The  variation  between  the  average 
densities  during  active  and  quiet  periods  is  a  factor  of  3 
and  related  to  the  decimetric  flux  (specifically,  the  10.  7- 
cm  radiation).  Jacchia' s  model  relates  csensity  directly, 
resulting  in  a  factor  of  3  for  all  altitudes.  Paetzold' s 
model  results  in  factors  greater  than  10  at  200  n  mi  that 
generally  increase  with  altitude.  In  addition  to  the  11-year 
cycle,  there  are  27- day  cycles  and  semiannual  and  annual 
cycles,  which  result  in  approximately  25-percent  varia¬ 
tions  at  100  n  mi  a«.>d  also  increase  with  altitude. 

c.  Magnetic  Storms.  These  are  generally  unpredictable,  but 
their  effects  are" relatively  short-lived,  lasting  for  only  a 
few  o*ys  The  effects  are  proportional  to  the  storm's 
intensity  and  can  vary  an  much  as  40- percent  at  100  n  mi 
and  much  more  at  higher  altitudes. 

Based  on  the  \>ove,  and  on  the  premise  that  100  n  mi  is 
the  principal  altitude  for  parking  orbits,  etc.  ,  tn~  nominal 
(mean)  atmospheric  density  model  was  conservatively 
chosen  as  the  1959  ARDC  (see  Table  E-l).  rt  car,  be 


o 


scaled  up  or  down  to  .include  the  average  solar  activity 
effects.  *  Combining  all  the  effects  of  denaity  and  drag 
coefficient  uncertainties,  no  as  to  estimate  a  standard 
deviation  for  purposes  oi  assessing  the  relative  impor¬ 
tance  of  aerodynamic  effects,  a  conservative  value  of 
0.2  (standard  deviation  of  EQOO)  is  recommended. 


Other  density  models  can  be  easily  input,  if  required,  for  which  Figure  E- 
presents  curves. 


( 


E- 1  "* 


-  iiaMc 


. —  -  AROC  1959 

- -  —  JACCHIA  I960 

F  *  ISO,  2  P.M.  LOCAL  TIME 
CONVERTED  FROM  TABLE  2,  REE  II 

- —  PAETZOLD  1962 

F»I50,  2  RM.  LOCAL  TIME 
USING  DATA  FROM  FIG.  3,  REF.  10 
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Figure  3.  Initial  Gyro  0-;entation 
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AS  SHOWN,  INITIAL  CONDITION  0ftlt»7*TI0N  OPTION  2 

? 
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Figure  5.  Coordinate  System  for  Initial  Condition  Errors 


DOWN  RANOE  IS  DEFINED  ALONG  THE  PROJECTION  OF  VELOCi.Y 
VECTOR  ONTO  HORIZONTAL  PLANE 

V,  IS  PITCH  COMPONENT  OF  VELOCITY  ERROR 

V?  IS  YAW  COMPONENT  OF  VELOCITY  ERROR 

&  IS  VELOCITY  MAGNITUDE  ERROR 

Figure  6  Coordinate  System  tor  Terminal  Condition  Errors 
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Table  G-l.  Error  Sources  . .  , . . .  G-3 

Table  G-2.  Orientation  and  Control  Data .  G-6 

Table  G-3.  Program  Constants .  G-9 


G-l 


*  These  definitions  are  con3istem  with  initial  condition  Option  2  and  platform  orientation 
Option  1. 


Table  G- 1 ,  Error  Sources  (Continued! 


* 


% 


i 


i 


i 


G-4 


Descriptive  of  single-  degree -of-freedom  gyro. 


Table  G-l.  Error  Sources  (Concluded) 


ee-of-freedom  gyro, 


Table  G-2.  Orientation  and  Control  Data 


Abort  time  for  reading  trajectory  tape 


Table  G-2.  Orientation  and  Control  Data  (Concluded} 
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Table  G-3.  Program  Constants  (Conversion  Factors) 


NOTE:  Underlined  numbers  only  are  program  constants 
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